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Abstract 


We calculate the independent helicity amplitudes in the decays B K* i+e- 
and B piu£ in the so-called Large-Energy-Effective-Theory (EEET). 
Taking into account the dominant 0{as) and SU{3) symmetry-breaking 
effects, we calculate various Dalitz distributions in these decays making use 
of the presently available data and decay form factors calculated in the QCD 
sum rule approach. Differential decay rates in the dilepton invariant mass and 
the Eorward-Backward asymmetry in B ^ are worked out. We also 

present the decay amplitudes in the transversity basis which has been used in 
the analysis of data on the resonant decay B K*J/Tp{-^ Measure¬ 

ments of the ratios Ri{s) = dTHiis){B K*i~^i~)/dTHiis){B pH-i'i), 
involving the helicity amplitudes Hi{s), i = 0,1, as precision tests of 
the standard model in semileptonic rare i?-decays are emphasized. We argue 
that i?o(s) and R-{s) can be used to determine the CKM ratio and 

search for new physics, where the later is illustrated by supersymmetry. 


Zusammenfassung 

Wir berechnen die unabhangigen Helizitatsamplituden der Zerfalle B 

und B piui in der sogenannten Earge-Energy-Effective-Theory 
(EEET). Unter Beriicksichtigung der dominierenden 0{as) und SU{3) 
Symmetrie-brechenden Effekte berechnen wir verschiedene Dalitz Distribu- 
tionen in diesen Zerfallen unter Einbeziehung der z. Zt. verfiigbaren Daten 
und Eormfaktoren, die mittels QCD Summenregeln berechnet wurden. Dif- 
ferentielle Zerfallsraten in der Dilepton-invarianten Masse und der Vorwarts- 
Riickwarts Asymmetric in B ^ wurden ausgearbeitet. AuBer- 

dem prasentieren wir die Zerfallsamplituden in der Transversalitatsbasis, die 
bei der Analyse der Daten des resonanten Zerfalls B K*J/ip{-^ 
benutzt wurde. Messungen der Verhaltnisse Ri{s) = dTHi{s){B 
K*/dVHi{s){B p^vt) mit den Helizitatsamplituden Hi(s),i = 

0, +1, — 1 werden als Prazisionstests des Standarsmodells in semileptonischen 
seltenen H-Zerfallen hervorgehoben. Wir diskutieren, daB Ro{s) und R-{s) 
benutzt werden kdnnen, um das CKM Verhaltnis |Eub|/|V)s| zu bestimmen, 
und um nach neuer Physik zu suchen. Eetzteres wird mittels Supersymmetric 
illustriert. 


1 




Contents 


1 Introduction 1 

2 Rare B Decays: Motivation and Methods 5 

2.1 The Flavour Sector in the Standard Model. 5 

2.1.1 Flavour changing neutral currents. 7 

2.2 The Effective Hamiltonian Theory. 9 

2.2.1 QCD improved corrections. 11 

2.3 6 —^ S 7 in the Effective Hamiltonian Theory . 14 

2.4 Eong-Distance Effects in Exclusive 5-decays. 17 

2.4.1 The Earge Energy Effective Theory (EEET). 18 

3 Exclusive B K*i+i- Decay in the SM 25 

3.1 Introduction. 25 

3.1.1 Kinematics. 27 

3.1.2 NEO-corrected amplitude for 6s . 27 

3.2 Hadronic matrix elements for 5. 32 

3.2.1 hard-spectator corrections. 35 

3.2.2 vertex corrections. 37 

3.2.3 form factors values. 39 

3.3 Decay Distributions in 5 —. 42 

3.3.1 Dalitz distributions. 44 

3.3.2 Dilepton mass spectrum and Eorward-backward asymmetry. 47 

3.3.3 Transversity Amplitudes for 5. 51 

3.4 Decay Distributions in 5 —> p . 57 

3.4.1 Dalitz distributions. 60 

3.4.2 Dilepton mass spectrum. 63 

3.5 Phenomenological Discussion on . 65 

3.5.1 Phenomenology. 67 

3.5.2 Model-independent analysis of Rb . 68 

3.6 Summary and Outlook. 71 

iii 




























4 Exclusive B K*i+i- Decay in SUSY 73 

4.1 Introduction. 73 

4.2 The Minimal Supersymmetric Standard Model . 74 

4.3 MSSM parameter space . 79 

4.4 Supersymmetric contribution to 6 —> s . 80 

4.5 Analysis in supersymmetry. 85 

4.5.1 SUGRA model. 85 

4.5.2 The MFV model. 87 

4.5.3 The Extended-MFV model. 89 

4.6 The ratios (s) and Ro{s) as probes of New Physics in B ^ 90 

4.7 Summary and Outlook. 92 

5 Summary & Future 94 

A Generalities 98 

A. 1 Input Parameters . 98 

A.2 QCD . 99 

A.3 Feynman Rules.100 

A. 4 Utilities.101 

B The large energy expansion 103 

B. 1 Feynman Rules in the Large Energy Limit of QCD .103 

B.2 The Eactorizabe corrections AFj.104 

B. 3 The Eunctions .105 

C B ^ K*i+i- in SUSY 108 

C. l The functions/i(a;).108 

C.2 The Auxiliary Eunctions ml, m|), 7712, m|) .109 

Bibliography .110 


IV 























Chapter 1 


Introduction 


Rare B decays have always played a crucial role in shaping the flavour structure of the Stan¬ 
dard Model [1] and particle physics in general. Since the first measurement of rare radiative 
B ^ K*'y decays by the CLEO collaboration [2] this area of particle physics has received 
much experimental [3] and theoretical [4] attention. In particular, flavour-changing neutral cur¬ 
rent (FCNC) 5-decays, involving the &-quark transition 6 —> (s, d) -f 7 and b —> (s, d) + 

(i = e, /i, r, u), provide a crucial testing grounds for the standard model at the quantum level, 
since such transitions are forbidden in the Born approximation. Hence, these rare 5-decays are 
characterized by their high sensitivity to New physics. 

In the standard model, the short distance contribution to rare 5-decays is dominated by the 
top quark, and long distance contributions by form factors. Precise measurements of these tran¬ 
sition will not only provide a good estimate of the top quark mass and the Cabibbo Kobayashi 
Maskawa (CKM) matrix elements [5] Vtd, Vts, Vtb, but also of the hadronic properties of 5- 
mesons, namely form factors which in turn would provide a good knowledge of the correspond¬ 
ing dynamics and more hint for the non-perturbative regime of QCD. 

Via the machinery of operator product expansion (OPE) and the renormalization group equa¬ 
tions (RGB) in the framework of an effective Hamiltonian formalism [6-10] (see section 2.2 for 
a discussion) 5e// ~ X] CiOi, one can factorize low energy weak processes in terms of pertur- 
batively short-distance Wilson coefficients Ci from the long-distance operator matrix elements 
< Oi >. The (new) vertices Oi, which are absent in the full Eagrangian, are obtained by in¬ 
tegrating out the heavy particles, namely the W and the top in the SM, from the full theory. 
Their effective coupling is given by the Ci, which characterize the short-distance dynamics of 
the underlying theory. 

The Wilson coefficient Gf®", reflecting the 6 —> 57 transitions, is actually well constrained* 
by the current precise measurement of the inclusive radiative 5 — Xg'y decays at the 5- 
factories. The current world average based on the improved measurements by the BABAR [11], 

*The modulus of the effective coefficient of the electromagnetic penguin operator in the SM agrees well with the 
experimental bounds, but there is no experimental information on the phase of C'y^^(m_B). 
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Decay Mode 

B Ke+e- 
B —> 

B Ki+i- 
B K*e+e- 
B —> 

B K*i+i- 
B XsC+e- 
B ^ Xsfi^fi~ 
B XJ+i- 


B (xlQ-^) 


BELLE [14] 


BABAR [18] 


± 0.06 

0.8+°;| ± 0.08 
0.58+°;}^ =b 0.06 


ri 9/1+0.49+0.14 
'^•^^-0.32-0.15 

l.SStlfs ± 0-26 

fi 7Q+0.24+0.11 
' °-0.20-0.18 


< 2.4 

< 1.2 
< 1.4 


1 70+0.87+0.48 
' O-0.72-0.49 

0.99+J;^^ ± 0.39 
1.68lg;| ±0.28 


5.0±2.3l};? 

7.9 ±2.1+2-0 

6.1 ± lAtli 


Table 1.1: Experimental results of the semileptonic rare B-decays [14,18]. 

CLEO [12], AEEPH [13] and BEEEE [14] collaborations, 

B{B ^ X, 7 ) = (3.43l°;f7) x 10-^ (1.1) 

is in good agreement with the estimates of the standard model (SM) [15-17], which we shall take 
as B{B —i- Xgy) = (3.50 ± 0.50) x 10“"^, and moreover, can exclude large parameter spaces of 
non-standard models. 

The b ^ s £+f“ (£+ = e+, /r+, r+) transition involves besides the electromagnetic penguin 
b sy* s also electroweak penguins b sZ^* —> s f+f“ and boxes. They give rise to 
two additional Wilson coefficients in the semileptonic B — s- Xsi^i~ decays, Cg and Ciq. 

A model independent fit of the short-distance coefficients Cg, Cg and Cio can be ob¬ 
tained, using the exclusive as well as the inclusive semileptonic (and radiative) rare 5-decays ex¬ 
perimental constraints on the corresponding branching ratio {^B (X^, K, K*)y ^, the various 

angular distributions and the Eorward-Backward (EB) asymmetry [19] in 5 —(X^, K, 
decays. They involve independent combinations of the Wilson coefficients, which allows the de¬ 
termination of sign and magnitude of Cf^^, Cg and Uio from data. On the other hand, these 
measurements are also of great help in studying that part of strong interaction physics which is 
least understood, the non-perturbative confinement interactions. 

Using the recent 5 —> X^y experimental result (1.1), with the new measurements of the 
semileptonic rare 5-decays recently reported by BEEEE [14] and BABAR [18] (see Table 1.1), 
it has been shown that the bounds on these Wilson coefficients are consistent with the SM, but 
considerable room for new physics effects are not excluded [20]. 

With the advent of the Eermilab booster BTeV (Eermilab) and EHCb (CERN) experiments 
at hadron colliders, and also the ongoing experiments at CEEO and the 5-factories, the semilep- 
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tonic rare decays of -B —(X^, K, will be precisely measured. On the theoretical side, 

partial results in next-to-next-to-leading logarithmic (NNLO) accuracy are now available in the 
inclusive decays B Xsi^i~ [21,22]. What concerns the exclusive decays, some theoretical 
progress in calculating their decay rates to NLO accuracy in the B {K*,p)'y [23-25], and 
to NNLO accuracy in B ^ [25] decays, including the leading Aqcd/^^b, has been 

reported. 

Making use of the exclusive semileptonic B ^ V (with V stands for a vector) theoret¬ 
ical improvements, obtained within the large-energy-effective-theory [26,27], we explore here a 
detailed phenomenological analysis of the exclusive B K* and B ^ p decays with 
i = e, p (since we neglected lepton masses in our calculation the result cannot be used in the r 
case) in the SM and supersymmetric theories . 

This thesis contains the following points [28,29]: 

• Using the effective Hamiltonian approach, and incorporating the perturbative improve¬ 
ments [25], we expressed the various helicity components in the decays B in 

the context of the Large-Energy-Effective-Theory. 

• As this framework does not predict the decay form factors, which have to be supplied from 
outside, we combined the B —> X *7 experimental constrains obtained in the NEO-EEET 

(K*) 

approach with the light cone QCD sum rule [30,31] to extract the EEET form factor 
and ^ in the large 7?^*-region s < 8 GeV^ (namely the EEET validity range). 

• We calculate a number of Dalitz distributions and the dilepton invariant mass distribution 

for the individual helicity amplitudes (and the sum) in B ^ We find that the 

NEO order corrections to latter distribution is significant in the low dilepton mass region 
(s < 2GeV2). 

• We show the 0{as) effects on the forward-backward asymmetry, confirming essentially 
the earlier work of Beneke, Feldmann and Seidel [25]. 

• We have compared the EEET-based transversity amplitudes in this basis with the data [32- 
35] currently available on B ^ K* J and find that the short-distance based 
transversity amplitudes are very similar to their long-distance counterparts. 

• Using the 5'(7(3)-breaking corrections, we relate the B p EEET form factors namely 

and with the B K* corresponding ones, and we implement the 0(Q;s)-improved 
analysis of the various helicity components in the decays B pine- We carry out in the 
context of the EEET a number of Dalitz distributions, the dilepton invariant mass distribu¬ 
tion for the individual helicity amplitudes (and the sum). 

• Combining the analysis of the decay modes B —> K*i~^i~ and B pine, we show that 

the ratios of differential decay rates involving definite helicity states, (s) and Ro{s), can 

be used for extracting the CKM matrix elements |14fe|/| V^^l. 
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• We investigate possible effects on these ratios from New Physics contributions, exemplified 
by representative values for the effective Wilson coefficients in the large-tan SUGRA 
models. 

Organization of the work 

An introduction to rare B decays and the methods used is given in Chapter 2, where we discuss 
the effective Hamiltonian theory, rare radiative B Xg'j decays and long-distance method in 
exclusive B decays. In Chapter 3 we investigate in details an helicity analysis of B ^ K* i+i- 
and B ^ p fvi in the SM. We present in the context of the NLO-LEET approach, various an¬ 
gular distributions and their uncertainties are worked out. Eurthermore, for the B 
decays, we project out the forward-backward asymmetry and the corresponding transversity ba¬ 
sis. Chapter 4 is devoted to the semileptonic rare B decay, by contrasting its antic¬ 

ipated phenomenological profile in some variants of supersymmetric models. After a review on 
the 6 —s decay in the MSSM, we propose to study the ratios Ro{s) and R-{s) as probes 
of new physics effects in B ^ using some generic SUSY effects. Einally, Chapter 5 

contains a summary and an outlook. Input parameters, Eeynman rules and utilities are collected 
in Appendix A. The large energy expansion with its machinery is presented in Appendix B. 
Appendix C contains various loop functions, introduced in SUSY. 




Chapter 2 


Rare B Decays: Motivation and Methods 


In this chapter we outline the flavour structure of the standard model (SM). We discuss the CKM 
mixing matrix and motivate the importance of studying flavour changing neutral current (FCNC) 
b ^ s transitions. We introduce the necessary tool to include QCD perturbative corrections in 
weak decays, the effective Hamiltonian theory. As an application of the former we discuss the 
6 —> S 7 decay as the most prominent example of a rare B decay. Finally we discuss in details 
the large energy quark expansion technique as the appropriate non-perturbative approach for the 
heavy-to-light transitions. 


2.1 The Flavour Sector in the Standard Model 


In the quark sector of the SM, there are six quarks organized in 3 families. The left-handed 
quarks are put into weak isospin SU{2)i doublets 



i=l,2,3 





( 2 . 1 ) 


and the corresponding right-handed fields transform as singlets under SU{2)l- Under the weak 
interaction an up-quark (with Qu = 2/3e) can decay into a down-quark (with Qd = — l/3e) and 
a boson. This charged current is given as 


jCC 


—— (m, c, t) 7/,Uckm 
y2sm.0w 


( 2 . 2 ) 


where the subscript L = (1 — 75)/2 denotes the left-handed projector and reflects the V — A 
structure of in the SM. Here the weak mixing (Weinberg-)angle 9w is a parameter of the 
SM, which is measured with high accuracy [36]. The so-called Cabibbo Kobayashi Maskawa 
(CKM) matrix Vckm [5] describes the mixing between different quark flavours. It contains the 
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angles describing the rotation between the eigen vectors of the weak interaction (g') and the mass 
eigen states (g) 




f d\ 

s' 

— Vckm 

s 

U'J 




(2.3) 


Symbolically, Vckm can be written as 


H:km 


\ Vtd Vts 


Vub \ 

Vcb 

Vtb / 


(2.4) 


In general all the entries are complex numbers, only restricted by unitarity VckmVckm = 
They are parameters of the SM and can only be obtained from an experiment. Note that only 
three independent real parameters and one phase are left after imposing the unitarity condition. 
Some parametrizations of Vckm can be seen in ref. [36]. 

A useful parametrization of the CKM matrix has been proposed by Wolfenstein [37] 




^Wolfenstein ~ 




A 

A 1 

AA ^(1 — p — i?7) — AA^ 




AA3(p 

AA 2 

1 


i7]) \ 




(2.5) 


The parameters A, A, p and the phase p are real numbers. A is related to the Cabibbo angle 
through A = sin6*c [36], which describes the quark mixing with 4 quark flavours. Since A ~ 
0.221, the relative sizes of the matrix elements in eq. (2.4) can be read off from eq. (2.5). As can 
be seen, the diagonal entries are close to unity and the more off-diagonal they are, the smaller 
is the value of their corresponding matrix elements. The parameter A has been determined from 
the decays b cii/i and B —> D*iu£, yielding A = 0.81 ± 0.04. The measurement of the ratio 
\Vub/Vcb\ = 0.08 ± 0.02 yields \/+ rf = 0.36 ± 0.09. Eikewise the mass difference AM^ = 
— M{B^p) ~ 0.50 (ps)“^ constrains the combination a/(1 — pY + p^. The observed 
CP-asymmetry parameter ex = 2.26 ■ 10“^ constrains p and p. The precise determination of the 
parameters p and p is a high and important goal, since it corresponds to two important questions: 


• Does CP hold in the SM ?? A non zero phase p 7 ^ 0 in the CKM matrix directly leads to 
CP violating effects. 


• The unitarity of the CKM matrix can be used to write down relations among its elements 
= ^ik-, i, /c = 1, 2, 3. There are 6 orthogonality equations possible (i 7 ^ 
k), and each can be represented graphically as a triangle, a so-called unitarity triangle 
(UT) [38]. The sides and angles of such an UT can be constrained by different types of 
experiments. Eor the UT given by the relation 


VudV;, + KW*: + v^bV:, = 0, 


(2.6) 
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there are 3 scenarios possible, which at present are not excluded experimentally and are a 
sign for new physics: 1. the UT does not close, i. e., 7^ 0’ where Oj denotes the 

three angles of the triangle. 2 . — 0 ’ the values of the are outside of their 

SM ranges determined by another type of experiment 3. ~ hut the values of 

the angles are inconsistent with the measured sides of the triangle. 

In the literature special unitarity triangles are discussed. A recent review over the present status 
on the CKM matrix and the unitarity triangle is given in [4]. 


2.1.1 Flavour changing neutral currents 

In the SM, the neutral currents mediated through the gauge bosons 7, g do not change flavour. 
Therefore, the so-called Flavour changing neutral currents (FCNC) do not appear at tree level and 
are described by loop effects. The subject of the present work is an analysis of such rare (FCNC 
mediated) B decays in the SM. The quarks are grouped into light {u, d, s) and heavy (c, b, t) 
ones in the sense, that the mass of a heavy quark is much larger than the typical scale of the 
strong interaction, Aqcd ~ 200 MeV. The sixth quark, the top, is too heavy to build bound 
states because it decays too fast. The special role of the 6 -quark is that it is the heaviest one 
building hadrons. We will not discuss the “double” heavy Be and concentrate on B = (bq) 
meson transitions with light q = u,d, s. Since the 6 -quark is heavy, the 5-system is well suited 
for a clean extraction of the underlying short-distance dynamics. Unlike the Ff-system, long¬ 
distance effects are expected to play a subdominant role in B decays except where such effects 
are present in a resonant form. 

The motivation to investigate b ^ s{d) transitions is to improve the knowledge of the CKM 
matrix elements and to study loop effects. For the latter the interest is large, since there is no tree 
level FCNC decay possible in the SM. The leading loops give the leading contribution and they 
are sensitive to the masses and other properties of the internal virtual particles like e. g. the top. 
They can be heavy and therefore can be studied in a rare B decay at energies which are much 
lower than the ones necessary for a direct production of such particles. The idea is to compare 
the SM based prediction for a rare B decay with an experiment. A possible deviation gives a hint 
not only for the existence, but also for the structure of the “new physics” beyond the SM. 

Further the 5-system can be used as a testing ground for QCD, to check perturbative and 
non-perturbative methods. One example is the decay 5 —^^ 7 , which can be described in the 
lowest order at parton level through 6 — 57 . As a 2-body decay, the photon energy in the 6 -quark 
rest frame is fixed: = (m^ — ml)/2mb for an on-shell 7. A possible non trivial spectrum 

can result from gluon bremsstrahlung 6 —syg and/or a non-perturbative mechanism, which is 
responsible for the motion of the 6 -quark inside the meson thus boosting the distribution. 

Some exclusive rare 5 decays have already been detected. The recent experimental observa¬ 
tions of the rare decay mode 5 — K*y have been determined by CLEO [32], and more recently 
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also by BABAR [34] and BEEEE [35] 


{ (4.55 ± 0.70 ± 0.34) x 10“^ [32], 
(4.39 ± 0.41 ± 0.27) x lO”® [34], 
(4.96 ± 0.67 ± 0.45) x lO'^ [35], 
and 

f (3.76 ± 0.86 ± 0.28) x lO”® [32], 
B{B+ K*+-i) = { 

1 (3.89 ± 0.93 ± 0.41) x lO'^ [35]. 


(2.7) 


( 2 . 8 ) 


However, the first observation of the rare 5-decay to the orbitally excited strange mesons has 
been reported by CEEO [32] and recently confirmed by BEEEE [14] with a branching fraction 
of 


B{B a:; (1430)7) 


(1.66lg;| ± 0.13) X 10-5 [32], 

(1.5;[|| ± 0.1) X 10-5 [14], 


(2.9) 


These important experimental measurements provides a crucial challenge to the theory. Many 
theoretical approaches have been employed to predict the exclusive B 77* (892)7 decay rate 
(for a review see [39] and references therein). On the other hand less attention has been devoted to 
rare radiative 5-decays to excited strange mesons [40-43]. Most of these theoretical approaches 
rely on non-relativistic quark models [40,41], HQET [42], relativistic model [43] and light cone 
QCD sum rules [44]. However there is a large spread between different results, due to a different 
treatment of the long distance effects. Thus, the difficulty with the exclusive mode is the large 
theoretical uncertainties due to the hadronic matrix elements, which has to be controlled. A large 
section is devoted to this issue in end of this chapter. 



✓ S 

/ \ 

/ \ 

I t 

_I_I 

h U, c, t s 


Eigure 2.1: A FCNC b ^ s diagram. 


A typical diagram for a virtual b ^ s transition is displayed in Eig. 2.1 from where the CKM 
couplings can be directly read off. The amplitude T is the sum over all internal up-quarks 

T= ^ ; A, = (2.10) 

i=u,c,t 

Using the CKM unitarity J2i=u ct^i = 0 the smallness of 14b yielding A„ \t, we arrive 
at 


T = 


At(Tb - T4 + A„(T„ - T4 ~ \t{Tt - TJ 


( 2 . 11 ) 
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for a 6 —> s amplitude in the SM. In the D-system the FCNC transition rates (c —> u) are mueh 
more suppressed due to an inbuilt GIM meehanism [45]. Here we have 

T = ^ V ■V*-T- 

i=d,s,b 

= V,tV:,(n-T,) + V^V;,{T,-T,), (2.12) 

in whieh the first term is CKM suppressed and the seeond one GIM suppressed sinee 

The SM rates in the eharm seetor for deeays sueh as 77, l'^l~ are out of 

reach for present experiments. If one nevertheless finds something in the rare charm sector, it 
would be a direct hint for the desired physics beyond the SM. 


2.2 The Effective Hamiltonian Theory 


As a weak decay under the presence of the strong interaction, rare B decays require special tech¬ 
niques, to be treated economically. The main tool to calculate such rare B decays is the effective 
Hamiltonian theory. It is a two step program, starting with an operator product expansion (OPE) 
and performing a renormalization group equation (RGE) analysis afterwards. The necessary 
machinery has been developed over the last years, see [6-10], [46] and references therein. 

The derivation starts as follows: If the kinematics of the decay are of the kind that the masses 
of the internal particles m* are much larger than the external momenta p, nif 3> p^, then the heavy 
particles can be integrated out. This concept takes a concrete form with the functional integral 
formalism. It means that the heavy particles are removed as dynamical degrees of freedom from 
the theory, hence their fields do not appear in the (effective) Eagrangian anymore. Their residual 
effect lies in the generated effective vertices. In this way an effective low energy theory can be 
constructed from a full theory like the SM. A well known example is the four-Eermi interaction, 
where the W -boson propagator is made local for -C (q denotes the momentum transfer 
through the W)\ 


9i^i' 




m 


1 q^ 
+ 


w 


m 


+ ... 


(2.13) 


w 


where the ellipses denote terms of higher order in 1/mw- * Performing an OPE for QCD and 
electroweak interactions, the effective Hamiltonian for a ECNC b ^ S'j transition in the SM can 
be obtained by integrating out W,t,4>. Up to 0(^^) it is given as: 


n.ffib ^ S7) = , (2.14) 

^ 2=1 

*We remark here that the original way was reversed; The main historical step was to extrapolate the observed 
low energy 4-Fermi theory in nuclear /3-decay to a dynamical theory of the weak interaction with heavy particle 
exchange. 
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where the weak couplings gw 


sin 9w 

Gf 

V2 

Gf = 


are collected in the Eermi constant 

9w 

1.16639 ■ 10“® GeV^ . 


Gf 


(2.15) 

(2.16) 


The on-shell operator basis is chosen to be [6,7] 

Oi = 4(sl«7^6i,„)(cl/37^cl/3), 

O 2 = 4 (sLa1fj.bLl3 ) (cLg'J^CLa ), 

O 3 = 4 (slo 7 ^ 6 l„) 

q=u,d^s,c,b 

O 4 = 4(slo7^6l/3) Y 

q=u,d,s,c,b 

O5 = 4 ( sl „ 7 ^& l „) Y 

q=u,d,s,c,b 

Oe = 4 ( sl „ 7 ^ 6 l / 3 ) Y i^RgY^Ra), 

q=u,d,s,c,b 

O 7 = a^^irribR + msYb 

Os = -^s^a^'^{mbR + msL)T^Y9G%: (2.17) 

where ttem = is the electromagnetic fine-structure constant, L{R) = 1/2(1 =f 75), = 

Ib/i, 7u\ and a, P are SU (3) colour indices. T“, a = 1... 8 are the generators of QCD, some 
of their identities can be seen in appendix A.2. Here G°'^’^ denote the electromagnetic 
and chromomagnetic field strength tensor, respectively. As can be seen from the operator basis, 
only degrees of freedom which are light compared to the heavy integrated out fields (kF, f, 0 ), 
remain in the theory. The basis given above contains four-quark operators Oi,,,q, which differ 
by colour and left-right structure. Among them, the current-current operators Oi and O 2 are 
the dominant four-Eermi operators. A typical diagram generating the so-called gluonic penguins 
Os ...6 is displayed in Eig. 2.2. The operators O 7 and Og are effective b ^ s'y, b ^ sg vertices, 
respectively. All operators have dimension 6 . Eor b —> sYi~ transitions the basis eq. (2.17) 
should be complemented by two additional operators containing dileptons. They are discussed 
together with their corresponding Wilson coefficients in chapter 3. 

The coupling strength of the introduced effective vertices Oi is given by the (c-numbers) 
Wilson coefficients Gpg)- Their values at a large scale /i = mw are obtained from a “matching” 
of the effective with the full theory. In the SM, the Gpmw) read as follows [47] 


Cl,3...6(mH/) 

= 0, 

(2.18) 

G2{rnw) 

= 1, 

(2.19) 

Gjimw) 

3a;^ — 2x^ —8x^ — + 7x 

4(a;-l)4 24(a;-l)3 ’ 

(2.20) 

Gs{mw) 

— 3a;^ —x^ -f 5a;^ -f 2x 

4(a;-l)4’^^+ 8(a;-l)3 ’ 

(2.21) 




11 


^ 


✓ \ 

/ \ 



Figure 2.2: A gluonic penguin diagram. 


with X = rn^/rn^. It is convenient to define ejfactive coefficients of the operators 

Oj and Os- They contain renormalization scheme dependent contributions from the four-quark 
operators (9i...6 in Tie// to the effective vertices inb ^ S'j and b —> sg, respectively. In the NDR 
scheme ^ , which will be used throughout this work, they are written as [7] 

Cf{p) = CM + QdC5{p) + QdN,CM , ( 2 . 22 ) 

Cfip) = Csip) + C,ip). (2.23) 

Here TV), denotes the number of colours, Nc = 3 for QCD. The above expressions can be found 



Figure 2.3: The diagram contributing to the one loop b ^ s'y, b ^ sg matrix element, respec¬ 
tively. 

from evaluating the diagram shown in Fig. 2.3. Contributions from (9i...4, which correspond to 
an ® like insertion, vanish for an on-shell photon, gluon, respectively. The Feynman 
rules consistent with these definitions are given in appendix A.3. 

2.2.1 QCD improved ctg corrections 

Our aim is now to include perturbative QCD corrections in the framework of the effective Hamil¬ 
tonian theory. This can be done by writing down the renormalization group equation for the 
t We recall that in the naive dimensional regularization (NDR) scheme the 75 matrix is total anti-commuting, i. 
e. {75,7 m} = 0, thus = 7^i?. 
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Wilson coefficients f 


(2.24) 


where 7 denotes the anomalous dimension matrix, i.e., in general the operators mix under renor¬ 
malization. Solving this equation yields the running of the couplings Cj(/r) under QCD from the 
large matching scale (here /i = mw) down to the low scale /i which is the relevant one 

for 6 -decays. Eq. (2.24) can be solved in perturbation theory = Aidas'. 


Iji 

a(/x) 


7® + ( 


167r2 


^ 7 ]!^ + • • 


IGtt^ 

r? 

9 ry /.a(1) 


IGtt^ 


Ci{fiy + 


(2.25) 

(2.26) 


The initial values of the above RGE are the Ci{mw), which in the lowest order in the SM are 
given in eq. (2.18-2.21). 

Eet us for the moment concentrate on the special case that 7 is a number. Then the lowest 
order solution is given by 


C(/i) = 

asijmw) 

V = -, 


(2.27) 

(2.28) 


which can be easily checked by substituting it into eq. (2.24). In the derivation we have used the 
QCD function, which describes the running of the strong coupling: 

Pig) = 9'~^9 = ~9iS 2^0 + + • • • , (2.29) 

d/i lOTT^ lOTT^ 

with its lowest order solution 


=i_ ( 230 ) 

aPfi) i + ^^f^ln(!<) 

We see that our obtained solution eq. (2.27) contains all powers of Q;s(/i) ln(^^). It is called lead¬ 
ing logarithmic (EEog) approximation and is an improvement of the conventional perturbation 
theory. In general such a QCD improved solution contains all large logarithms like n = 0,1,... 
(here with g, = mb) 


a-(m,)ln-(^), 

mw 


(2.31) 


where m = n corresponds to EEog. A calculation including the next to lowest order terms is 
called next to leading order (NEO) and would result in a summation of all terms with m = n — 1 


With Ci = Ci{g, g) we have equivalently 
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C^{|x) 

/i = mw 

/i = 10 GeV 

fj, = 5 GeV 

/i = 2.5 GeV 

C'l 

0 

-0.161 

-0.240 

-0.347 

C's 

1 

1.064 

1.103 

1.161 

6*3 

0 

0.007 

0.011 

0.016 

6*4 

0 

-0.017 

-0.025 

-0.035 

a 

0 

0.005 

0.007 

0.010 


0 

-0.019 

-0.030 

-0.046 


-0.196 

-0.277 

-0.311 

-0.353 

^eff 

-0.098 

-0.134 

-0.148 

-0.164 


Table 2.1: Leading order Wilson coefficients in the Standard Model as a function of the renor¬ 
malization scale jj,. 


and so on. In the following we use the 2-loop expression for as^jj) which can be always cast into 
the form 


(/^) 


47r 

/3o 




(2.32) 


With Nf = 5 active flavours (note that we integrated out the top) and SU{Nc = 3) the values of 
the coefficients of the /9 function are 


Po — 


23 

y ’ 



(2.33) 


They are given in appendix A.2 for arbitrary Nc and Nf. The strong scale parameter Aqcd = 
A^qcd^'^ is chosen to reproduce the measured value of Q:s(/i) at the pole. 

We recall that in LLog the calculation of the anomalous dimension and the matching con¬ 
ditions at lowest order, Cf°^(mvi/) is required. In NLO a further evaluation of higher order 
diagrams yielding cl^\mw) is necessary and in addition the hadronic matrix elements (Oi) 
have also to be known in 0{as). 

In a general theory and also in the one relevant for rare radiative b decays given in eq. (2.14), 
the operators mix and the matrix 7 has to diagonalized. In the latter case the (8 x 8 ) matrix 7 *^°^ 
has been obtained by [8,9] and the running of the Cj(/i) in LLog approximation cannot be given 
analytically. The LLog solution for the Wilson coefficients ready for numerical analysis can be 
taken from [48]. We display the Ci for different values of the scale jj, in Table 2.1. As can be 
seen, there is a strong dependence on the renormalization scale fx, especially for Ci and 
Other sources of uncertainty in the short-distance coefficients Ci are the top mass and the value 
of as{mz)- We keep them fixed to their central values given in appendix A.l. 

Here a comment about power counting in our effective theory is in order: As can be seen from 
Fig. 2.3 with an external photon, the insertion of four-Fermi operators generates a contribution to 
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b —> s'f, which is also called a “penguin”. It is a 1-loop diagram, but unlike “normal” perturbation 
theory, of order q;°. To get the a] contribution, one has to perform already 2 loops and so on. 
That means, the calculation of the EO(NEO) anomalous dimension matrix was a 2(3)-loop task. 

A comprehensive discussion of weak decays beyond leading logarithms can be seen in ref. [46]. 
The main results of the NEO calculation in B ^ Xg'j decay will be given in section 2.3. 

The advantages of the effective theory compared to the full theory can be summarized as 
follows: 

• The effective theory is the more appropriate way to include QCD corrections. Earge loga¬ 
rithms like \n{fx/mw) multiplied by powers of the strong coupling as{fJ,), which spoil the 
perturbation series in the full theory, can be resummed with the help of the RGE. 

• On the level of a generic amplitude A = (Tie//) ~ the problem can 

be factorized into two parts: The short-distance (SD) information, which can be calcu¬ 
lated perturbatively, is encoded in the Ci, and it is independent of the external states, i.e. 
quarks or hadrons. The long-distance (ED) contribution lies in the hadronic matrix ele¬ 
ments. Both are separated by the renormalization scale /r. Of course the complete physical 
answer should not depend on the scale /i, truncating the perturbation series causes such a 
remaining dependence, which can be reduced only after including higher order terms or a 
full resummation of the theory. 

• As long as the basis is complete, the effective Hamiltonian theory enables one to write 
down a model independent analysis in terms of the SD coefficients Ci. This is true for SM 
near extensions like the two Higgs doublet model (2HDM) and the minimal supersym¬ 
metric model (MSSM). Here one can try to fit the Ci from the data [49]. However, new 
physics scenarios like, e.g., the left-right symmetric model (ERM) require an extended op¬ 
erator set [50-52]. Wilson coefficients in the 2HDM and in supersymmetry (SUSY) can 
be seen in ref. [53] and ref. [54], respectively. 

2.3 6 ^ §7 in the Effective Hamiltonian Theory 

The effective Hamiltonian theory displayed in the previous section is applied to 6 —> sy tran¬ 
sitions. Several groups have worked on the completion of the EEog calculation [8,9]. The 
anomalous dimension matrix at leading order 7 *^°^ and the lowest order matching conditions 
(eq. (2.18-2.21)) govern the running of the EEog Wilson coefficients, denoted in this and only 
this section by cf'* (/r), to separate them from the NEO coefficients. We discuss the improvement 
of the theory in B ^ obtained from NEO analysis. In the remainder of this work we treat 
the Wilson coefficients Uj, i = 1,... 8 in EEog approximation. 

In the spectator model, the inclusive B branching ratio in EEog approximation can 
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be written as 


B{B 


Xsj) = Bsi 


T{b^ 57) 
r(6 — ceVe) 


|-^t| SOem 

I Kb 12 7r/(mc) 


|^(o)eff 


(/^)P 


(2.34) 


where a normalization to the semileptonic decay B Xciiy^ to reduce the uncertainty in the 
6 -quark mass has been performed. Here Bsi denotes the measured semileptonic branching ratio 
and the phase space factor^ firric) for T{B ^ Xciue). 

As the branching ratio for B Xg'y is mainly driven by several effects can be 

deduced: 


• Including LLog QCD corrections enhance the branching ratio for B about a factor 

2 — 3, as can be seen in Table 2.1 (here denoted by Ci{jj,)) and changing the scale from 
/i = mw down to /i ~ nib. 

• While the sign of is fixed within the SM, i.e. negative, it can be plus or minus in 

possible extensions of the SM. A measurement of B{B Xs'j) alone is not sufficient to 
determine the sign of or in general, the sign of Cfresulting from possible higher 

order calculations. 

• The strong scale dependence of the value of causes serious problems in the 

accuracy of the LLog prediction. Varying the scale between ^ < /i < 2mb , results in an 
error in the branching ratio of ±25% [7,55]. 


Because of the last point the NLO calculation was required. Several steps have been necessary 
for a complete NLO analysis. Let us illustrate how the individual pieces look like: At NLO, the 
matrix element for b ^ sj renormalized around fx = nib can be written as [7]: 


M.{b ^ 57 ) 



\tD{0’j{nib)) tree 1 


(2.35) 


with 


D = c7±) + 


asirub ) 

An 


Y. 




nr 


+ Cr%Ar„ 

fl 


(2.36) 


The ViT are computed in ref. [57]. They contain the b s'yg bremsstrahlung corrections [6,56] 
and virtual corrections to the Oj matrix element [57]. Especially the latter with an O 2 operator 
insertion demands an involved 2-loop calculation, see Figs. 1-4 in [57], where the corresponding 
diagrams are shown. It is consistent to keep the pieces in the parentheses in eq. (2.36), which are 
multiplied by as{mb), in LLog approximation. 

Now Cf^^(/r) has be be known at NLO precision, 

. (2.37) 

_ 47r 

§For the semileptnic B Xcih'i decay, the phase space factor read as: f{mc) = 1 — 8mc^ + Srh^ — rh^ — 

24toc^ Irirric (toc = me/mb) . 
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As this job consists out of two parts, the work has been done by two groups: The 0{ag) anoma¬ 
lous dimension matrix was obtained in ref. [58], which required the calculation of the residue 
of a large number of 3-loop diagrams, describing the mixing between the four-Eermi operators 
and Oj^s- The second part, the NEO matching at /i = mw has been done in ref. [59] and 
confirmed in ref. [60]. The NEO calculation reduces the /i = 0{mb) scale uncertainty in varying 
fjt, in the range ^ < /r < 2mb drastically to ±4.3% [61] and suggests for B a scale 

/r = ^ as an “effective” NEO calculation through 

T{B ^ X, 7 )io(/i = ^) ~ r(5 ^ • (2.38) 

As a final remark on scale uncertainties it should be noted that in the foregoing the top quark 
and the W have been integrated out at the same scale /r = mw, which is an approximation to be 
tested. It is justified by the fact that the difference between as{mw) and as{mt) is much smaller 
than the one between as{mw) and as{mb) The authors of [61] analysed the dependencies 
on both the W matching scale fiw = and the one at which the running top mass is 

defined: mtifit) and fiw % fJ't- Similar to the mb scale they allowed for nw, j^t a possible range: 
^ < /U-x < 2ma: where x = W,t. Their findings are that the /i^y, /r* uncertainty is much smaller 
(namely ±1.1%, ±0.4% at /i ~ m;, in NEO, respectively) than the uncertainty in the scale around 
mb and therefore negligible. 

Concerning the exclusive b ^ S'j transitions, the situation is more complicated. Eor a generic 
radiative decay B F 7 , where F = P,V,S,A,T and Ta stands for pseudoscalar, vector, 
scalar, axial-vector, tensor and pseudo-tensor respectively, one defines the corresponding exclu¬ 
sive branching ratio as following: 

B{B^P{S)^) = 

B{B^V{A)-^) = 

X 

B{B^T{Ta)i) = 

X 

where ms {mp) and tb are the 5-meson mass (the generic F-meson) and life time respectively, 
whereas Pf (g^) is the so-called transition form factor, which will be given in section 3.2.3. 

A good quantity to test the model dependence of the form factors for the exclusive decay is 
the ratio of the exclusive-to-inclusive radiative decay branching ratio: 

^ B{B ^ 1/(A) 7 ) 

_ B{B X, 7) 

^Using eq. (2.32) and the input parameters in Table A.l, we have as{mw) = 0.12, as('tnt) = 0.11 and 

Ois{mb) = 0 . 21 . 


Ctp 


Gl 


tb- 


\\b?ml\cf{mb)?F^^^\Qf 


327r4 




rriB^ 

Mml\cf{mbrF^^^-\0f 

5 


m 


B 


m; 


'T{Ta) 


1 ± 


m. 


PTa) 


m 


B 


(2.39) 


(2.40) 


(2.41) 
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and 




(l-^y(A)) (} + ^ViA)) 

(1 - ml/rrib'^)^ (1 + m^jm^) ’ 


'^T(Ta) 


B{B ^ T{Ta) 7) 

B{B ^ X, 7) 

^r('rA)(o)2 j 

^^t{Ta) (1 - (1 + mljmb^') ’ 


(2.42) 


(2.43) 


where rhp = mpImB- With this normalization, one eliminates the uneertainties from the CKM 
faetor \ and the short distanee Wilson eoeffieient Cf^^{mb). Thus, we are left in eqs.(2.42) 
and (2.43) with unknown form faetors Ff'{0), whieh have to be eomputed using some non- 
perturbative methods, whieh will be presented in the fortheoming seetion. 


2.4 Long-Distance Effects in Exclusive B-decays 


After having witnessed in explieit terms the short-distanee (eoeffieients) of the OPE, we will now 
turn to the long-distanee (operator matrix elements) eontributions in exclusive B Deeays. 

In general, 5-mesons transitions ean be measured inelusively over the hadronie final state 
or exelusively by tagging a partieular light hadron (typieally a Kaon for B K*^* transition). 
The inclusive measurement is experimentally more diffieult but theoretieally simpler to interpret, 
sinee the deeay rate is well and systematieally approximated by the ealeulation of quark level 
proeesses. However, the theoretieal diffieuly with exelusive deeay modes is usually due to their 
nonperturbative nature eneoded in their hadronie form faetors. 

For a 5-meson deeay into a pseudosealar meson (5), the eorresponding form faetors are 
defined by the following Lorentz deeomposition of bilinear quark eurrent matrix elements: 


{P{p)\q'lMB{PB)) = f+{q^ 


{Pb + p)t, - 


m\ — rnj> 


% 


+ /o(g^ 


— nip 


Qt., (2.44) 


{P{p)\q(T^^q’'b\B{pB)) = [q‘^{pB + p)^^ - {ml - ml) q^] , (2.45) 

tub nip 

where is the 5 meson mass, mp the mass of the pseudosealar meson and q^ = Pp ~ P^ is 
the four-momentum transfer . The relevant form faetors for 5 deeays into veetor meson (K) are 
defined as 


{V{p, e*)\q-ff,b\B{pB)) 
{V{p, €*)\qj^j5b\B{pB)) 


2V{q^ 


rriB + mv 


P^PBi 


2inivAQ{q^ 


■ q 


qp + i {mp + mv) Ai{q‘^) 


(2.46) 

e* ■ q 


q" 
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(I/(p, e*)\qa^uq''h\B{pB)) 
{V (p, £*) \qa^^,^5q''h\B{pB)) 


-i A-2{q^ 


■ q 


rriB + mv 


(Pb + p)^, - 


m\ — Triy 


% 


2iTi{q^) e^^p^e*"' p'^bP'', 

T2{q^) [(m| - my) - [e* ■ q) {pb + p)p\ 


+ T3(g2) {e*-q) 




m\ — my 


{Pb+p)^ 


(2.47) 

(2.48) 

(2.49) 


where my {s) is the mass (polarisation vector) of the vector meson, Aq, 1,2(7^) and Ti_2,3(7^) are 
defined respectively as the semileptonic and the penguin form factors. 

Clearly, in order to compute these form factors ^/+,o,r(7^), Aq, 1,2(7^) and Ti^2,3(7^)j one 
is forced to use some theoretical methods such as the Heavy Quark Effective Theory (HQET), 
the Earge Energy Effective Theory (EEET), QCD sum rules, Eattice QCD or Quark Models. 
Needless to say, all these non-perturbative methods have some limitations. Consequently the 
dominant theoretical uncertainties in the exclusive modes reside in these form factors. 

Among all these theoretical approaches, it has been shown recently, that an adequate tool 
to describe heavy-to-light i?-transitions is the so-called Earge Energy Effective Theory (EEET). 
Since we are dealing in this analysis with B {K*, p) transition, we will focus in the following 
on the EEET approach, more appropriate for our work. 


2.4.1 The Large Energy Effective Theory (EEET) 

Although the HQET [62] has permitted a great succes in the description of heavy-to heavy 
semileptonic decays such as B it fails unfortunately in its description of the heavy- 

to-light decays, such as B ^ M transitions, where M stands for a light mesonIL 

The EEET was first introduced by Dugan and Grinstein [26] to study factorization of non- 
leptonic matrix elements in decays like B D^*'>7r, D(*) p..., where the light meson is emitted 
by the kE-boson. Eater on, Charles et al. [27] have established this formalism for semileptonic 
and radiative rare 5-decays, such as B ^ Tiiui, B and B —> K*'j. They have 

shown that to leading order all the weak current P P{V) matrix elements can be expressed in 
terms of only three universal form factors. However, the EEET symmetries are broken by QCD 
interactions and the leading 0(0;*) corrections in perturbation theory are known [25,63]. 


Heavy-light form factors at large recoil 

Eet us switch to the system under consideration, the 5-mesons. Since the 6-quark (inside the 5- 
meson) is heavy, i .e. mi, ~ 4.6 GeV >> Aqcd ~ 200MeV, it will transmit all its momentum 
to the light quark (inside the final light meson P{V)) with a large energy Eptyp in almost the 
whole physical phase space except the vicinity of the zero-recoil point: 


Ep{y) — 


mp 

2 


1 - 


rn 


m 


+ 


P{V) \ 


B 


m 


(2.50) 


B 


II a bound states of light quarks: u, d and s. 
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(c) (d) 

Figure 2.4: Dijferent contributions to the B —» PiV) transition, where u stands for light quark 
{u, d, s). (a) Soft contribution (soft interactions with the spectator antiquark q' are not drawn), 
(b) Hard vertex renormalisation. (c,d) Hard spectator interaction. 

This assumption holds in the limit where sueh transitions are dominated by soft gluon exehange, 
i. e. the 6-quark and the light one must interaet with the speetator quark (and other soft degrees 
of freedom) exelusively via soft exehange, as it is shown in Fig. 2.4a. 

To work out the large-reeoil symmetry eonstraints on the soft form faetor, one uses a teeh- 
nique familiar from HQET [64-67], in writing the heavy-to-light eurrent [gF^lgcD = [QrJ'bv]eS- 
The form faetors at large reeoil ean be ealeulated within the following set up of the LEFT: 

• The heavy 6-quark momentum is written as pb = rribV + k, where A; is a small residual 
momentum of order Aqcd and v denotes the veloeity of the meson with momentum pb = 
tubv whieh at rest isv = (1, 0, 0, 0). 

• The heavy 6-quark is then deseribed by its large spinor field Q{x) eomponent by{x) = 

• The light g-quark momentum earries SLPq = Ep{y)n- + k' momentum, where k' ~ Aqcd 
and the light-like veetor n_ (nf = 0) is parallel to the four-momentum p = Ep{y)n- of 
the light meson P{V). 

• The light g-quark is deseribed by the large eomponents g„(a;) = q{x) of its 

quark spinor field q{x). Here rip = 2v — n_ is another light-like veetor with n+ ■ n_ = 2 
and Eg ^ Ep(y) is the energy of the light quark. 

Following the above deseription, the form faetors at large reeoil are then represented by (defining 
L = P,V) 


{L{ELn-)\qJPK \B{mBv)) = ti [A l{El) Ml^ Mb] , 


(2.51) 
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where 

-Mp(v) = -25 

-Mb = ^(- 75 ), (2.52) 

with e the polarisation vector of the vector meson. The function Al{E) contains the long¬ 
distance dynamics, but it is independent on the Dirac structure T of the current, because the 
effective lagrangians (see Appendix B.l) does not contain a Dirac matrix. The most general 
form Al{E) can take is therefore 

Al{El) = aiiiEL) + a2L{EL) j + a^iiEL) i- + a^iiEL) i-j-, (2.53) 

but the projectors AIl, -Mb imply that not all the aniEi) are independent. Accounting for these 
projectors, the most general form is 

Ap{Ep) = 2Epi^^\Ep), (2.54) 

Ay{Ev) = Eyi_\i^P{Ey)-ti^P{Ey)Y (2.55) 

with a conveniently chosen overall normalisation. It follows that the three pseudoscalar meson 
form factors are all related to a single function (Ep) and the seven vector meson form factors 
are all related to two unknown functions, ^^^\Ey) and ,^|^|^^(i?y). The latter two functions are 
chosen such that only ^^^\Ey) contributes the form factors for a transversely polarised vector 
meson and only ,^|^|'^^(i?y) contributes the production of a longitudinally polarised vector meson. 
Performing the trace in Eq. (2.51), we obtain 

{P{p)\qYb\B{pB)) = 2Ep^<^^\q^)n^_, 

{P{p)\qcT^''qub\B{pB)) = 2iEpi^P\q^) [[rriB - Ep) n^_ - (2.56) 

for pseudoscalar mesons, and 

{V{p,e*)\qYh\B) = t2Ey {q^) , 

{V{p,e*) \qj'"j5b\ B{pb)) = 2Ey |^l'^^(g^) [e*^ - (e* ■ v)n^_] + iP{q^)^{e* • u)n((| , 

{V{p,e*)\qa^''qub\B{pB)) = 2EymB iP P) elvpU-a, 

(E(p, £*) B{pb)) = -i2EyiP{q^) - e*^np 

-t2Ey^P{q^)p^{e* ■ v) - n^_v^) , (2.57) 

for vector mesons, in agreement with Ref. [27]. Comparing Eqs. (2.44)-(2.49) with Eqs. (2.56)- 
(2.57), we find the following form factor relations: 


fAf) = ^''■>( 9 "). 


( 2 . 58 ) 
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Symmetry-breaking corrections to the LEET form factors 

We have just seen the LEET effect in describing the exclusive heavy-to light semileptonic decays 
by reducing the number of independent form factors from ten to three. However theses symme¬ 
tries are broken by factorizable and non-factorizable QCD corrections, worked out by Beneke et 
al [25,63], 

While the form factors obtained in Eqs.(2.58)-(2.67) are a straighforward evaluation of the 
soft contributions in Eig. 2.4a, their (9(a^)-corrections are originated from the two following 
processes: 

• the vertex corrections (Eig. 2.4b) 

• the hard scattering corrections (Eig. 2.4c-d) 

The vertex corrections are a straightforward calculations using standard techniques, in contrast 
to the hard scattering ones where one makes use of the two-particle light-cone distribution am¬ 
plitudes of the B meson and the light meson (more details can be found in Ref. [63]). 

Einally, having these C>(Q;s)-corrections at hand, the form factors defined in Eqs.(2.58)-(2.67) 
get modified as follows [63]: 


( 2 . 68 ) 




22 


CHAPTER 2. RARE B DECAYS: MOTIVATION AND METHODS 


/o(g^) 


2Ei 


niB 


-{''“'tf) T + ^ |2-2y]l.+ 


An 




tub 


as 

An 


1.4 


as Cb 
A n 

+ 2y 


A/o, 


+4^ap, 

An 


(2.69) 

(2.70) 


for the form factors of pseudoscalar mesons and 
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for the form factors of vector mesons. The abbreviation y stands for 


3^ = 


2E 


P(V) 


In 


2E 


P{V) 


mp — 2Ep(^y) mp 


(2.71) 

(2.72) 

(2.73) 
AA 2 , (2.74) 

(2.75) 

(2.76) 

(2.77) 

(2.78) 


Moreover in Eqs. (2.68)-(2.77), the form factors receive a further additive correction from the 
interaction with the spectator quark, indicated by AEi (and can be found in Appendix B.2). Its 
general form reads as 


AEi ^^B®Ti®^L 

where Ti is a hard-scattering kernel convoluted with the light-cone distribution amplitudes of the 
B meson and the light meson L. Thus, we can summarize the 0 (q;s)-LEET corrections by the 
following, tentative, factorization formula for a heavy-light form factor at large recoil, and at 
leading order in 1/mp: 

Mq^) = + (2.79) 

where ^^{Ep) is the soft part of the form factor, to which the EEET symmetries discussed above 
apply and 4 = 1-1- 0{as) is the hard vertex corrections. At this stage, we have seen that the 
0 (q;5 )-EEET corrections can be absorbed into the redefinition of the corresponding EEET form 
factor 4 
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Figure 2.5: Various next-to-leading order eontributions to the B ^ V^* matrix elements. 

However, concerning the B —> V"-transition, there exist further corrections at order o;*, orig¬ 
inating from four quark operators and the chromomagnetic dipole operator in the weak effective 
hamiltonian, which cannot be expressed in terms of form factors, i.e. matrix elements of the 
type {V\sTb\B), and will be presented in section 3.2. Sample Feynman diagrams are shown in 
Fig. 2.5e-g, compared to the diagrams in Fig. 2.5a-d, which do assume the structure of form 
factor matrix elements and it will be discussed in section 3.2. 


Soft-collinear contributions to the LEET form factors 

An important unresolved question in strong interaction physics concerns the parameterization of 
power-suppressed long-distance effects to hard processes that do not admit an operator product 
expansion (OPE). For a large class of processes of this type the principal difficulty arises from 
the presence of collinear modes, i.e. highly energetic, but nearly massless particles. 

Recently Bauer et al. [68] have claimed that the missing collinear gluons in the LEET do 
not allowed this effective theory to reproduce the Infrared (IR) physics of QCD [69]. Thus an 
effective theory is able to reproduce correctly the infrared physics of QCD at one loop, only by 
including both collinear and soft gluons [70]. 

Happily, Bauer et al. [68] and Beneke et al. [71] have formulated separately the heavy-to- 
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n+p ~ 1, ~ 


Eigure 2.6: Kinematics of heavy quark decay into a single cluster of collinear and ultrasoft 
particles. 


light soft-collinear effective theory by taking into account the collinear and ultrasoft particles 
(see Eigure 2.6), missing in the EEET approach. They found independently that the presence of 
collinear gluons does not spoil the relations among the soft form factors, therefore establishing 
the corresponding results in the large energy limit of QCD [63]. 

Since the numerical effect of the collinear gluons contributions are negligeable [71] on the 
EEET form factors (defined in section 2.4.1), we will not consider them in our work. 




Chapter 3 


Exclusive B K*t^l Decay in the SM 


This chapter contains a comprehensive helicity analysis of the B and the B piV 

decays in the so-called Large-Energy-Effective-Theory (EEET). Taking into account the domi¬ 
nant 0{as) and SU{3) symmetry-breaking effects, we calculate various double and single dis¬ 
tributions in these decays making use of the presently available data and decay form factors 
calculated in the QCD sum rule approach. As precision tests of the standard model in semilep- 
tonic rare 5-decays, we propose a model independent extraction of the CKM matrix elements 
Kb\/\Vts\. 


3.1 Introduction 

Elavour changing neutral current (ECNC) decays b ^ s 'j and b ^ s are governed in the 
SM by loop effects. They provide a sensitive probe of the flavour sector in the SM and search 
for physics beyond the SM. In the context of rare B decays the radiative mode b ^ s 'j has been 
extensively discussed in chapter 2. 

In this chapter we address the exclusive semileptonic B and the B piV de¬ 
cays with in the EEET framework. Since we are neglecting finite lepton masses 

we cannot apply our results to the r-case. The theoretical study of the exclusive rare decays 
proceeds in two steps. Eirst, the effective Hamiltonian for such transitions is derived by cal¬ 
culating the leading and next-to-leading loop diagrams in the SM and by using the operator 
product expansion and renormalization group techniques (for a review see [72] and references 
therein). Second, one needs to evaluate the matrix elements of the effective Hamiltonian between 
hadronic states. This part of the calculation is model dependent since it involves nonperturba- 
tive QCD. Many theoretical approaches have been employed to predict the exclusive radiative 
B — {K, K*) t^t~ {i = e, p, t) decays. Most of them rely on QCD sum rules [73,31], quark 
model [74], lattice-constrained dispersion quark model in [75] and perturbative QCD [76]. Re¬ 
cently, an updated analysis of these decays has been done in [20] by including explicit 0{as) 
and Aqcd/^^6 corrections. Concerning, the lepton polarizations in the B K*{-^ Ktt) 
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decay in terms of the helicity amplitudes without the explicit 0{as) corrections, was 

undertaken in a number of papers [77-83]. In particular, Kim et al. [79,80] emphasized the role 
of the azimuthal angle distribution as a precision test of the SM. Eollowing closely the earlier 
analyses, we now calculate the 0{as) corrections in the LEET framework. 

Concentrating on the decay B the main theoretical tool is the factorization 

Ansatz which enables one to relate the form factors in full QCD (called in the literature A* (g^) {i = 
0,1, 2), U(g^), Tj(g^) {i = 1, 2, 3)) and the two LEET form factors and.^ll^^(g^) [25,63]; 

fkiq^) = ® Tfc ® , (3.1) 

where the quantities Ci {i =±, ||) encode the perturbative improvements of the factorized part 


Ci = cf^ + —cf) + ..., 

TT 


and Tfc is the hard spectator kernel (regulated so as to be free of the end-point singularities), 
representing the non-factorizable perturbative corrections, with the direct product understood as 
a convolution of Tk with the light-cone distribution amplitudes of the B meson ($b) and the 
vector meson (<l>y). With this Ansatz, it is a straightforward exercise to implement the 0{as)- 
improvements in the various helicity amplitudes. The non-perturbative information is encoded 
in the LEET-form factors, which are a priori unknown, and the various parameters which enter 
in the description of the non-factorizing hard spectator contribution, which we shall discuss at 
some length. The normalization of the LEET form factor ^ (g^) at g^ = 0 is determined by the 
B —> 76*7 decay rate; the other form factor \q‘^) has to be modeled entirely for which we 
use the light cone QCD sum rules. This input, which for sure is model-dependent, is being used 
to illustrate the various distributions and should be replaced as more precise data on the decay 
B —piug becomes available, which then can be used directly to determine the form factors 
^( 7 ^) and ,^ 11 ^ ^ (g^), taking into account the SU (3)-breaking effects. 

This chapter is divided roughly into three parts. The first one (this section up to and including 
section 3.3), contains an introduction to B ^ decay, basic definitions and the 0{as) 

improvements to the b ^ s matrix elements in the LEET framework. It is mainly devoted to 
the analysis of the double and single angular distributions for the individual helicity amplitudes, 
and their sum, and the Eorward-Backward (EB) asymmetry. In doing that we have shown the 
systematic improvement in 0{as) and 1/E in the exclusive radiative B decay, using 

the large energy expansion (LEET). Eurther, we carry out in the so-called transversity basis, the 
LEET-based transversity amplitudes (both in the LO and NLO accuracy), and compare them to 
the currently available data. 

The second part discussed in section 3.4, describes a helicity distributions analysis of the ex¬ 
clusive semileptonic B decay in the LEET. We display the various helicity components, 

Dalitz distributions, and the dilepton (yet) invariant mass, making explicit the 0{as) corrections. 
The estimates of the B p LEET form factors ^^^(s) and ^||^^(<s), which are scaled from their 
B ^ K* counterparts incorporating SU(3)-breaking, are also displayed here. 
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Finally the last part deals with subseetion 3.5.2, is devoted to the determination of the ratio of 
the CKM matrix elements | Vub | /1 ft® | from the ratio of the dilepton mass speetra inB ^ piu^ and 
B —> deeays involving definite helieity states. In partieular, we show the dependenee of 

t+t- Ids 

the ratio RAs) = —- (i = 0,-1) involving the helieity-0 and -1 eomponents, on the 

^ Ids 

CKM matrix elements | Vuh I /1 V'ts |. 

3.1.1 Kinematics 

We start with the definition of the kinematies of the exelusive semileptonie B V ■, t I'f) 
deeays, 

B{pb) V{pv,ev) +i{p+) + {i,iye){p-) , (3.3) 

where the index V stands for the eorresponding veetor meson. We define the momentum transfer 
to the lepton pair and the invariant mass of the dilepton system, respeetively, as 

(3.4) 

(3.5) 

The dimensionless variables with a hat denotes normalization in terms of the 5-meson mass, 
e.g., 

s = \ = {P++P-f, 
my 

mv = -, 

tub 
TTT-e 

m = —, 

rriB 

ete., where my (rrii) is the eorresponding veetor meson mass (lepton mass). Sinee we are dealing 
in our analysis with the two lepton generations, namely = e^, /i^, one ean negleet their finite 
masses. Thus, the sealed variable s in the 5 —C ug) deeays is bounded as follows, 

(m^ ~ 0) < s < (1 — myY ■ (3.9) 


(3.6) 

(3.7) 

(3.8) 


q = + p_ , 

s = . 


3.1.2 NLO-corrected amplitude for h ^ s 

Next, the explieit expressions for the matrix element and (partial) branehing ratios in the de¬ 
eays b ^ s are presented in terms of the Wilson eoeffieients of the effeetive Hamiltonian 
obtained by integrating out the top quark and the kF^ bosons. 


Heffib 


■) = Heffib ^ S + 7) - -^KsVtb [CMOg + CioOio] , 


(3.10) 




28 


CHAPTER 3. EXCLUSIVE B K*i+i- DECAY IN THE SM 



cT 

C'2 





LL 

-0.257 

1.112 

0.012 

-0.026 

0.008 

-0.033 

NLL 

-0.151 

1.059 

0.012 

-0.034 

0.010 

-0.040 



'-'8 

C^g 

C'lo 

/^NNLL 

^NNLL 

'-"10 

LL 

-0.314 

-0.149 

2.007 

0 



NLL 

-0.308 

-0.169 

4.154 

-4.261 

4.214 

-4.312 


Table 3.1: Wilson coefficients at the scale /i = 4.6 GeV in leading-logarithmic (LL) 
and next-to-leading-logarithmic order (NLL) [25]. 

where Heff{b —s + 7) together with the operators O 1...8 and their corresponding Wilson coeffi¬ 
cients Cfjj) [7,6] can be seen in section 2.2. The two additional operators involving the dileptons 
Og and Oio are defined as: 

TT 

ry _ 

Oio = (3-11) 

TT 

A usual, CKM unitarity has been used in factoring out the product VfgVth- The Wilson coef¬ 
ficients are given in the literature (see, for example, [9,48]). They depend, in general, on the 
renormalization scale /i, except for Ciq. With the help of the effective Hamiltonian in eq. (3.10) 
the matrix element for the decay b si+i- can be factorized into a leptonic and a hadronic part 
as, 

M{b ^ si+n = V*VtJCg [s^^Lb] + Cio [s7^T6] 

V 27r t 

-2rhbCf siaf^^jRb | . (3.12) 

where we neglect the rris mass. Here and in the remainder of this work we shall denote by 
mb = rribiif) the MS mass evaluated at a scale /r, and by mb^poie the pole mass of the 6-quark. 
To next-to-leading order the pole and MS masses are related by 

mbiti) = mb,pole (^1 + 31n^-4 -f 0(a^)^ . (3.13) 

Since we are including the next-to-leading corrections into our analysis, we will take the Wilson 
coefficients in next-to-leading-logarithmic order (NLL) given in Table 3.1. 

The K* meson subsequently decays to K and tt, with effective Hamiltonian 

"^e// ~ 9K*Ktt{PK — Pn) ' ^K*, (3.14) 

with gK*Kit denotes the strong coupling of Tf-mesons, to P-wave pion. In the following analysis, 
we neglect the masses of leptons, kaon and pion. Then the final 4-body decay amplitude can be 
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written as the sum of two amplitudes, 


= Mr + Ml , 


( 3 . 15 ) 


where A4r and Ml denote respectively the left and right helicity amplitudes in the dilepton 
system; and they can be written in a compact form, 

Mr = —^VtbVt*sgK*Kn - {a-Rdfii/ “ ^Rp^q^ + icR e^yapP°‘q^) (3.16) 

V 2 TTS 


(,PK Pn^a-) 


— m|-* + irriK*'^ K* 

Ml = ^VtbV^lgK*K7T(iLYh) {aLg^iu - bLP^Qu + icL e^uagp'^q^) (3.17) 
V 2 TTS 


gvcx _ 


'(.PK Pn')ai 


p2 - 77i|.^ + iniK*^K* 

with p = pr + p.^, and the auxiliary functions aR,bR, cr and aL,bL,CL can be expressed as 

i(mB + rriK* 


<^L/R 

bi/R 


2 rrib mB\/s i 
i 


s mB(±Cio - C9)Ai(s) + 4Ti(s) mb{mK* - mB)EK* 
4 Ti(s) rub (-m| + m^.) Er* 


rUb tub (rn% ~ ^k*) ^ 


+mBs{ - 2 mb< Ti(s) + %(s) 


(3.18) 

(3.19) 


mB 

2Er* 


T2(s) } + 2)2(5) (±Cio - C9)(mij - tur *) 


cl/r 


mb (mB + mR*) ^/s 


2 Ti(s) mb (mB + mR*) + (=fCio + Cg) 5(^(5) 


(3.20) 


Note that our conventions for aL,R, bL,R and cl,r are slightly differents from those defined by 
Kim et al. in ref. [79] by a factor of 1/y/s. The form factors V, Aq, Ai, A 2 , Ti, T 2 and T 3 have 
already been introduced in section 2.4 and their numerical values will be presented below. The 
functions Ti, 7^ and 7^ are related to the so-called penguin form factors, and will be defined in 
the next section. 

With the help of the above expressions, the differential decay width becomes. 


_ fr _ 

dp‘^ ds d cos 9 r d cos 9 + dcj) 


2^/X 

128 X 2567 r 6 m| 


(\Mr\^ + \Ml?). 


(3.21) 


with 


A 



m 


K* 


m 


K* 



(3.22) 


Here we introduced the various angles as: 9 r is the polar angle of the K meson momentum in 
the rest system of the K* meson with respect to the helicity axis, i.e. the outgoing direction of 
K*. Similarly 6 ^+ is the polar angle of the positron in the dilepton rest system with respect to the 
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helicity axis of the dilepton. And 0 is the azimuthal angle between the planes of the two decays 
K* ^ Kn and 7 * —> And then, 


\M 


i?l 


^F-rr Tr* (^errJFlb 


an? { {P ■ <lf - {P ■ Q? - 


(p2 - + {'mK*TK*y 

{s - q'^)p'‘^ 


+ 2Re{aRb*ji) {-{p ■ qfip ■ q) + (p ■ q)ip ■ q')ip ■ g)} 

+ l&fiP I {p ■ qf{p' ■ qf - ip ■ qTip' ■ qY - ^ J p^{p> ■ qY 


+ \cr\ { -{q'pqp'){q'pqp') - 


s — q‘ 


,/2 


-{pqp') ■ (pqp') 


+ 2Im{bRC*ji) {p ■ q) {p ■ q) {q'pqp') + 2Im{cRa*^) {p ■ q) {q'pqp') 
- 2Im{bRa*jf){p ■ q){q'pqp') - 2Re{cRa*j^){qp'q') ■ (p'pq) 

+ 2Re{bRC%){p' ■ q)(^) ■ (Fm)] , 


(3.23) 


and 


\M 


L - 


^F-rr Tr* (^errJFlb 

—fF Vtb Vts9K*KTT - 

V2 ns 


(p2 - + {mK*TK*y 


a^r \ (P ■ qY - (P ■ qY - 


(s — g'^)s 


+ 2Re{aLbl) {-{p ■ qYip ■ q) + {p ■ q){,P ■ q){p' ■ g)} 

+ WI{p> ■ qYip>' ■ qY - (p ■ q'Yip' ■ qY - ^ J pYp ■ qY 


+ Yl\ S -{q'pqp'){q'pqp') - 


s — q 


/2 


-{pqp') ■ {pqp') 


2Im{bLcl){p ■ q'){p' ■ q){q'pqp') + 2Im{cLal){p' ■ q){q'pqp') 
2Im{bLal){p' ■ q){q'pqp') + 2Re{cLal){qp'q') ■ {p'pq) 
2Re{bLcl){p' ■ q)(^) ■ {Vm)] , 


(3.24) 


where* {ABC)^ = e^a/s-yA’^B^C'^, (ABCD) = Eaq^ysA^B/^C^D^, p' = px — P-k, and q' = 

p+ — P-. Comparing with we see that the signs of the corresponding last three 

terms are opposite to each other. We can simplify the expression by introducing the helicity 
amplitudes, defined as. 


H, 


-L 

^(±,0) 

(±1,0) - 

tR* 

-R 

A±fl) 

(±1,0) - 

Sr* 


{aLObiu - biPb^qu + icR p'^q'^) , 

{aRQtiu - bRP^q^ + icR p'^q^) , 


(3.25) 


where the auxiliary functions aip, bip, cl,r are given in Eqs. (3.18)- (3.20), and we define the 


“We use Tr(7“7^7'>'7‘^75) = +4ze“^''''^. 
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following polarization vectors: 


^K* — 

: 0, 

1, 

h 

0 

)/V2, 

Cr* = 

: 0, 

1, 

-h 

0 

)/V2, 

4* = 

^ ms ’ 

0, 

0, 


)/niK* 

4 = 

: 0, 

1, 

-h 

0 

)/V2. 

s = 

: 0, 

1, 

+i, 

0 

)/V2, 

6° = 

7 

y ms ’ 

0, 

0, 

+ -s 

V 

)/75. 


Substituting them into Eq. (3.25), we obtain the following helicity amplitudes, 


Hifis) 


{(^L,R + 

{0'L,R — Cl,R"/\) 

P-Q , ^L,R^ 

—a-L,R - 1 = H- -p 

mK*Vs mK*VS 


(3.26) 


(3.27) 


where p.q = —s)/2. From now on, we will omit for simplicity the index 1 in the helic¬ 

ity amplitudes^ in Eqs. (3.27). Using these equations, we can get the results for Eqs. (3.23,3.24), 
whose sum makes the decay angular distribution oi B ^ K*{-^ K7r)i^i~, 


dp^ ds dcosOx dcos6j^ dcj) 64 x 8 ( 27 r)®m|[(p 2 — 

X {4cos^6'i^ sin^6'+(|i7o^(s)|^-f |i7o^(s)n 

+ sin2 0^(l + cos^ e+){\H^{s)\^ + \H^{s)\^ + \HY{s)\^ + \HY{s)\Y 

-2siiYeKsiiYe+ [cos2(j)Re{H^{s)HY*{s) + 

-sin20 Im{H^{s)HY*{s) + H^{s)H!i*{s))] 

- sin 26K sin 26+ [cos(j)Re{H^{s)H^*{s) + HY{s)H^*{s) + H^{s)H/^*{s) + 

- sinci) ImiH^{s)H^*{s) - HYis)H^*is) + H^is)H/^*is) - H^is)H/^*is))] 

-2 sin^ 6 k cos 6+{\H^{s)\‘^ - \H^{s)\‘^ - \H+{s)\‘^ + \H^{s)\Y 

+2sin0+ sin20i^ [cos0i?e(i7:[^(s)77o^*(s) - HY{s)H^*{s) - H^{s)H/^*{s) + H^{s)H/^*{s)) 

- sin0 ImiH^{s)HY*is) + HYis)H^*is) - H^is)H/^*is) - H^is)H/^*is))] } • (3.28) 


From the decay angular distribution presented above, it turns out that the main theoretical 
difficulty in evaluating this quantity is the estimate of non-perturbative part located in the helicity 
amplitudes (see Eq .3.27). Henceforth, we will see in the next section our estimate of the related 
hadronic matrix elements. 


^We will refer to hYYY) respectively as H^^is) and hY^{s). 
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3.2 Hadronic matrix elements for B 


In this section, we present our estimates of the non-perturbative effects on the exclusive 
B —decay, which are described by the matrix elements of the quark operators in 
Eq. (3.12) over meson states, and can be parameterized in terms of form factors. 

Eor the vector meson K* with polarization vector the semileptonic form factors of the 
{V — A) current are defined as 


{K*{p,e*)\{V - A)f,\B{pB)) = -i e* {ruB + mK*)Ai{s) +i {pb + p)f,{e*PB) —^ 7 ^^^— 

ruB + rriK* 

+i qti{e*PB) ("^ 3 ( 5 ) - Ao(s)) + e*''p%p'^ —^ ^ 9 ^ 

Note the exact relations: 


Asis) 

Ao{0) 

{K*\dpAf^\B) 


mB + rriK* . . , 
-^i(«) ■ 

Ztuk* 

213 ( 0 ), 

2 mx*(e*Pi?)Ao(s). 


mB - rriK* 
2mK* 


^2(’S), 


(3.30) 


The second relation in (3.30) ensures that there is no kinematical singularity in the matrix element 
at s = 0. The decay B K*Ei~ is described by the above semileptonic form factors and the 
following penguin form factors: 

{K*{p,e*)\Cfs(Tp^q'"{l + 'y5)b\B{pB)) = i ep^p^e*''p%p'^ 2Ti{s) (3.31) 

+ ^2(s)|e*(m| -m^*) - (e^ps) { pb + v )).^ 

+ B{s){e*pB) - 2 1^2 iPB+P)A • 

The matrix element decomposition is defined such that the leading order contribution from the 
electromagnetic dipole operator Oj reads 7)(s) = Cf^Tj(s) + ..., where Ti{s) denote the tensor 
form factors. Including also the four-quark operators (but neglecting for the moment annihilation 
contributions), the leading logarithmic expressions are [84] 


Ti(s) = CfT,{s) + Y{s) 


Us) = CfT2{s) + Y{s) 


Us) = CfUs) + Y{s) 


2mb{mB + mx*] 
s 

2mb{mB — rriK*, 


V{s), 


mB - mx* , / ^ mB + mx* 


2mb 


^lof-s) — 


2mh 


Ai(s) 


(3.32) 

(3.33) 

(3.34) 


with Cf = Ct- Us/S - dQ/g - 20C'5/3 - SOUe/g = Ct- (4^3 - C' 5 )/g - (4^4 - C'6)/3, 
and the function Y (s) represents the one-loop matrix element of the four-Eermi operators [48,9], 
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see Fig. 3.1. It is written as: 


Y (s) — h{s, TTic) (3C*i + C*2 3“ 3C*3 3“ C*4 3- 3C*5 + Coj 

~ — h{s, mb) (4 (C *3 + C 4 ) + 3C*5 + Cq) — — h{s, 0) {C^ + 3 C 4 ) 

3- - ^ 2^3 3- 2(3*4 3- "^^5^ ’ (3.35) 


where the “barred” coefficients Ci ( for i=l,..., 6 ) are defined as certain linear combinations of 
the Ci, such that the Ci coincide at leading logarithmic order with the Wilson coefficients in the 
standard basis [85]. Following Ref. [25], they are expressed as : 


C. = ic. 

C2 = C'2-^C'i, 

0 

C 3 = <3*3 — - (3*4 + 16 (3*5 — - (3*6, 

6 3 

^4 = ^C'4 + 8C'6, 

^5 = C'3-^C'4 + 4C'5-^C'6, 

6 3 

^6 = ^6*4+ 2 a. (3.36) 

The functions 


h{s,mq) 


and 

h{s, 0 ) 



(2 3- z) \f\; 


arctan ■ 


1 < 


Vz — 1 


In 


1 3“ \f\ — z 


it: 

y 


Z > 1 , 
(3.37) 
z < 1 , 


8 4, s 4 . 

^ - 9 'V + 9 


(3.38) 


are related to the basic fermion loop. (Here is defined as Am^/s.) Y (s) is given in the NDR 
scheme with anticommuting 75 and with respect to the operator basis of [15]. As can be seen 
from the above equations, internal 6 -quarks ~ h{s,mb), c-quarks ~ h{s,mc) and light quarks 
q, (with rrig = 0 for g = u, d, s) ~ h{s, 0 ) contribute to the function Y(s); only the charm loop 
involves the dominant “current-current” operators Oi and O 2 . 

Since Cg is basis-dependent starting from next-to-leading logarithmic order, the terms not 
proportional to h{s, mg) differ from those given in [85]. The contributions from the four-quark 
operators Oi-q are usually combined with the coefficient Cg into an “effective” (basis- and 
scheme-independent) Wilson coefficient 


Cf(s) = C9 + F(5). 


(3.39) 
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The effective Wilson coefficient C 9 ^^(s) receives contributions from various pieces especially 
from the cc states^. However the contribution given below is just the perturbative part. 



Eigure 3.1: The Feynman diagram responsible for the four-Fermi-operator contribution (depicted 
by the blob) to the operator Og. 


We have seen in the previous chapter that in the Large Energy Effective Theory framework, 
one can relate the seven B K* form factors to only two universal quantities [27,63], namely 
^ and \ Adopting this formalism, the various form factors appearing in (3.32)-(3.34) 

simplify to 


T,is) ^ T^is) = f 

Us) 

Us) 


{K*)i 


Cf Si + 


2mbmB 


Y{s) 


2Ek* 

ms 

ms 


2E 


Us), 

ns) s T,(s) ^-rp\s) 


K* 


C 7 U 2 +^Y{s)Ss 


2mh 


(3.40) 

(3.41) 

(3.42) 


where Ek* refers to the energy of the final state 7f*-meson (see Eq. (2.50) in section 2.4.1). The 
factors Si are defined such that 5* = 1 + 0{as). The 0(Q;s)-corrections represent the next-to- 
leading terms related to these form factors in the LEET and they can seen from Eqs. (2.73)-(2.77). 

At next-to-leading order, the invariant amplitudes 71 ||(s), which refer to the decay into a 
transversely and longitudinally polarized vector meson (virtual photon), get contributions both 
from factorizable corrections as well from the non-factorizable ones, and they read respectively 
[25] 




47r 
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fBU,± 


^1 = eil 


K 

rn 


K 


E 


niB 


C',® + 
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± 

O-sCp 


doj 

^ ^B,±U / du^K*,U)T±,±{,u,uj), 


UJ 


An 


C, 


( 1 ) 


^This effect will not be treated here since the LEET symmetry is restricted to the kinematic region in which the 
energy of the final state meson scales with the heavy quark mass. Eor B K*£'^£~ decay this region is identified 


as s ~ 8 GeV . 
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Figure 3.2: Leading contributions to ^*\Hes\B). The circled cross marks the 
possible insertions of the virtual photon line. 

j du^K*,\\{u)T\\^±{u,u) . (3.43) 

Here Cp = 4/3, Nc = 3, = 1, S|| = mK*/E]^, fK*,\\ denotes the usual K* deeay eon- 

stant fx* and /x*,± refers to the (seale-dependent) transverse deeay constant defined by the 
matrix element of the tensor current. The leading-ordercoefficient Ca'^ follows by comparison 
with Eqs. (3.40) and (3.42) setting 6i = 1. The term Ta,±{u,u}) and (a =T, ||) represent 
respectively the hard-scattering and the form factor corrections which will be discussed below. 


3.2.1 hard-spectator corrections 


The hard scattering functions ±(m, u) (a =_L, 11) in Eq. (3.43) is expanded as : 


Ta,±{u,u) 



agCp 

47r 



(3.44) 


To compute the leading-order term we have to compute the weak annihilation amplitude of 
Eigure (3.2-c), which has no analogue in the inclusive decay and generates the hard-scattering 
term u) in (3.43). To compute this term we perform the projection of the amplitude on 

the B meson and K* meson distribution amplitude as explained in [63]. The four diagrams in 
Eigure (3.2-c) contribute at different powers in the 1/mh expansion. It turns out that the leading 
contribution comes from the single diagram with the photon emitted from the spectator quark 
in the B meson, because this allows the quark propagator to be off-shell by an amount of order 
the off-shellness being of order ml for the other three diagrams. Hence the result of 
the leading-order term u) reads [25] 




mBCO AmB 
mBOJ — s — le mb 


(C'3 + 3C'4). 


(3.45) 

(3.46) 


The hard scattering functions in (3.44) contain a factorizable term from expressing the full 

( K*) 

QCD form factors in terms of , related to the -correction to the 6i in Eqs. (3.40), (3.42) 
above. We write -f The factorizable correction reads [63]: 



2m B 
uEx* ’ 


(3.47) 
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Eigure 3.3: Non-factorizable contributions to K*\Hqq\B). The circled cross marks the possi¬ 
ble insertions of the virtual photon line. Diagrams that follow from (c) and (e) by symmetry are 
not shown. Upper line: hard spectator scattering. Lower line: diagrams involving a B ^ K* 
form factor (the spectator quark line is not drawn for these diagrams). 



Cf 


+ 


s 

2mbmB 


Y{s) 


= 0 


2m\ 


(3.48) 

(3.49) 


The non-factorizable correction is obtained by computing matrix elements of four-quark opera¬ 
tors and the chromomagnetic dipole operator represented by diagrams (a) and (b) in Eigure (3.3), 
using the projection on the meson distribution amplitudes. They read as [63]: 






Acd Cl** tub 


eut±{u, me) (C 2 -f C 4 - Cq) 


u + us/m\ 2mb 
+ ed t± {u, mb) (C 3 -f C 4 - Ce - Amb/ ms C 5 ) -f Cd t± {u, 0 ) C 3 


e„t|| {u, me) (C2 + C4 - Ce) + Cd ty {u, mb) (C3 -f C4 - Ce) 


0 

mB 
mb 

+ edt\\{u, 0) C3 
mBOO 


mBOJ — s — ie \-u + us/m 


8 C, 


eff 


B 


6m 


B 


mb 


(3.50) 

(3.51) 

(3.52) 


h{um\ -f us, me) (C 2 + C 4 -f Cq) + h{um% -f us, mb) (C 3 -f C 4 -f Cq) 


+ h{um\ -\- us, 0 ) (C 3 -f 3 C 4 -f 3Cq) — — (C 3 — C 5 — ISCg) 


(3.53) 


Here C|** = Cg + C 3 - C 4/6 + 2 OC 5 - lOCg/S = Cg + (dCg - Cq)/3, e„ = 2/3 (e, = -1/3), 
Cq is the electric charge of the spectator quark in the B meson and h{s, mq) has been defined 
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above. The funetions ta{u,mg) arise from the two diagrams of Figure (3.3-b) in whieh the 
photon attaehes to the internal quark loop. They are given by [63]: 


t±{u,mq) 




hifUq) + (^Bo{umB + us, niq) - Bo{s, m^) j, 

2^5 r / \ um% + us f ^ 2 \ n / \3 

Ji(mg)+ - [Bo{umq^ + us,mg) -Bo{s,mg)), 

uEk* u \ / 


(3.54) 

(3.55) 


where Bq and Ji are defined as 

Boi^XjUig) = —2 a/ Arus/x — 1 arctan 
2m, 


and 


Ji(m„) = 1 + 


x± = -± 


uirrin — s) 


m. 


^/Ams/x — l’ 

Li(a;+) + Li(a;_) - Li{y+) - Li{y_) 


4 um% + us 


1 , 

, y± = - ±\ - -^ 1 


Li(a;) = In -—^ ln(l — x) — — + Li 2 

X 6 


2 

a; 


4 s y ’ 


x — 1 


(3.56) 

(3.57) 


(3.58) 

(3.59) 


The eorreet imaginary parts are obtained by interpreting as m^ — ie. Closer inspeetion shows 
that eontrary to appearanee none of the hard-seattering funetions is more singular than 1 /u as 
M 1. It follows that the eonvolution integrals with the kaon light-eone distribution are eonver- 
gent at the endpoints. 

The limit s —0 (E^* ms/^) of the transverse amplitude is relevant to the deeay B 

K*-y. The eorresponding limiting funetion reads 


t±{u,mq)\ 


_ 4 

s=0 — = 
U 



2m2 

um% 


Li{x+) + Li{x-) 



(3.60) 


In the same limit the longitudinal amplitude develops a logarithmie singularity, whieh is of no 
eonsequenee, beeause the longitudinal eontribution to the B K* i+i- deeay rate is suppressed 
by a power of s relative to the transverse eontribution in this limit. 


3.2.2 vertex corrections 

The next-to-leading order eoeffieients in (3.43) eontain a faetorizable term from expressing 

the full QCD form faetors in terms of \ related to the a^-eorreetion to the 5i in Eqs. (3.40), 

(3.42). In writing = Ca^ + C^a^\ the faetorizable eorreetion reads [25]: 

Cf = Cf (^4 In ^ - 4 - , (3.61) 

C® = -Cf (^4 In ^ - 6 + 4L^ + ^ F(s) (2 - 2l) 

" \ /i^ / 2mb V / 


(3.62) 
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b Oi_2 s b Oi_2 s b Oi_2 s 



Eigure 3.4: Complete list of two-loop Feynman diagrams for b —> sy* associated with the op¬ 
erators Oi and O 2 . The fermions (b, s and c quarks) are represented by solid lines; the curly 
lines represent gluons. The circle-crosses denote the possible locations for emission of a virtual 
photon. 


with 


L = 


mt 


In I 1 — 


mt 


(3.63) 


Note that the brackets multiplying Cf® include the term 3 ln(m^//i^) — 4 from expressing the 
MS quark mass in the definition of the operator Of in terms of the b quark pole mass according to 
(3.13). The non-factorizable correction is obtained by computing matrix elements of four-quark 
operators and the chromomagnetic dipole operator represented by diagrams (c) through (e) in 
Eigure (3.3). 

The matrix elements of four-quark operators require the calculation of two-loop diagrams 
with several different mass scales. The result for the current-current operators Oi 2 is presented in 
[22] as a double expansion ins/ml and me/mb. Since we are only interested in small s, this result 
is adequate for our purposes. Eor that note that only the result corresponding to Eigure (3.4a-e) 
of is needed for this calculation. The 2-loop matrix elements of penguin operators have not yet 
been computed and will hence be neglected. Due to the small Wilson coefficients of the penguin 
operators, this should be a very good approximation. The matrix element of the chromomagnetic 
dipole operator [Eigure (3.3-c)] is also given in [22] in expanded form. The exact result is given 
in Appendix B.3. All this combined, we obtain 


CfC 


(nf) 


m ^eff p(7) 

- (^ 2^2 - Og i^g 


2mbmB 


C2F^^^ + 2Ci 




(9) 


CfC, 


(nf) 


, ruB 


eff Tn{7) 


+ CTF, 


2mb 


C2F!f^ + 2^1 


/ i?(9) 




Fr + gFf ) +C|"F, 


(9) 


eff 17(9) 


FCTF, 


■yeff 17(9) 


(3.64) 


(3.65) 


The quantities F^f''’ are given in Appendix B.3, or can also be extracted from [22] 

in expanded form. In expressing the result in terms of the coefficients Ui, 2 , we have made use 
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A 2 

^0 

V 

Ti 

T 2 

Ts 

m 

0.294 

0.246 

0.412 

0.399 

0.334 

0.334 

0.234 

Cl 

0.656 

1.237 

1.543 

1.537 

1.575 

0.562 

1.230 

C2 

0.456 

0.822 

0.954 

1.123 

1.140 

0.481 

1.089 


Table 3.2: Input values for the parameterization (3.69) of the B K* form factors. Renormal¬ 

ization scale for the penguin form factors Ti is p, = mb [31 ]. 

of + F 2 ^ /Q = 0. We also substituted Cs by taking into account a subset of penguin 
contributions. 

3.2.3 form factors values 

In the description of exclusive 5-decays hadronic matrix elements < X\Oi\B > are involved^. 
However, in order for these quantities to become available, it is necessary to confront the fact 
these hadrons are color bound state objects. While understood in principle, the non-perturbative 
nature of these bound states makes problematic the extraction of precision information about the 
exclusive 5-physics. To explore them one faces a daunting theoretical challenge to evaluate first 
the corresponding form factors. 

This is not a problem which has been solved in its entirely, nor is it likely ever to be. Rather, 
what is available is a variety of theoretical approaches and techniques, appropriate to a variety 
of specific problems and with varying levels of reliability. While approaches which are based 
directly on QCD, and which allow for quantitative error estimates, are clearly to be preferred, 
more model-dependent methods are often all that are available and thus have an important role 
to play as well. 

Concerning our guess on the corresponding B ^ K* form factors (see Eqs. (3.29) and (3.31)) 
we have combined roughly two theoretical approaches to compute them: 

• First we have used the LEFT symmetry to reduce the number of independent form factors 
from seven to two universal ones, namely \s) and ^(s) (see Eqs. (2.61)-(2.67)). 

• At large recoil, namely s = 0, the normalization of the FEET form factor ^(0) is 
determined using the 5 —> 5^*7 experimental constraint on the corresponding NEO-FEET 
estimates. Thus, the magnetic moment form factors turns out to be in the range [23-25] 

T;f^*)(0) = 0.28 ±0.04. (3.66) 

^where X is any meson (with mass mx) 
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(K*) 

Thus, from Eqs. (3.66) and (2.65), the numerical value for the large recoil momen¬ 

tum is defined as 


^(_^*)(0) = 0.28 ±0.04. (3.67) 

• The second universal LEET form factor ^(s) has to be modeled entirely from some 
approximate methods. Eor that we have used a non-perturbative approach, the so-called 
Light-cone sum-rule approach^ [86,87], based on the approximate conformal invariance 
of QCD. While in principle this technique is rigorous, it suffers in its current practical 
implementations from a degree of uncontrolled model dependence (for a review see [90] 
and reference therein). Using the result of [30] for B K* form factors, presented 
in Table 3.2, which include NLO radiative corrections and higher twist corrections up to 
twist four. The result turns out to be: 

^^*(0) = 0.31 ±0.04. (3.68) 


• To extrapolate these form factors, namely \ at small recoil (large values of s), we use 
the following parametrization (with its coefficients listed in Table 3.2): 


F{s) = F(0) exp(cis ± C 2 S^). 


(3.69) 


Using the ingredients described above, we show in Eig. 3.5 the corresponding LEET form factors 
{s) and .^11 \s). Note that the values used by Beneke el al. in ref. [25] are very different 

of the ones used by us. This descripency is related to the fact that their choice is based on the 
QCD sum rules estimates for ^ = 0.35 ± 0.07 and ^ = 0.49 ± 0.09. On the other hand, 
our values are somewhat lower than the corresponding estimates in the lattice-QCD framework, 
yielding [91] Tf*(0) = 0.32^°;°^, and in the light cone QCD sum rule approach, which give 
typically T^* (0) = 0.38 ± 0.05 [31]. (Earlier lattice-QCD results onB ^ 76*7 form factors are 
reviewed in [92].). 

Einally, we have to keep in mind that such a descripency reflect after all our poor knowledge 
of this part of QCD, namely the non-perturbative QCD. Consequently, we can anticipate the fact 
that the long distance uncertainty in our analysis will be the dominant one. 


^The method of light-cone sum-rules was first suggested for the study of weak baryon decays in [86] and later 
extended to heavy-meson decays in [87]. It combines the traditional QCD sum rule method [88] with the twist 
expansion characteristic of hard exclusive processes in QCD [89]. 
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Figure 3.5: LEET form factors \s) for B —> K* l~. The two columns denoted by [AS] 

and [BES] represent, respectively, our ^(s) and the ones used by Beneke et al. in ref [25]. 
The central values are represented by the dashed curves, while the bands reflect the uncertainties 
on the form factors [28]. 
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3.3 Decay Distributions in B ^ 

Using all the machinery presented until this point we are able now to do the corresponding 
analysis for specific transitions. In this section, we present our Helicty analysis of the B 
decay. In the LEET Limit, the helicity amplitudes (3.27) are expressed as: 




- 4 Ti{s) rub {niB - rriK*) {ms + Ek* 


2 mb ms (niB + rnK*)^/s L 
+ (±Cio - Cg) me {ruB + mK*Y s Ai(s) 
+2mBA/A|2Ti(s) mb {mB + mx*) + (=fCio + Cg) sV{s 


(3.70) 


hB^(s) 


4 Ti(s) mb {mB - nix*) {mB + mx*f Ex* 


2 mb mB {mB + mx*)y/s 
+ (±Cio - Cg)mB{mB + mx*Y s Ai(s) 
-2mB'/\\‘2Ti{s) mb{mB + mx*) + (=fCio + C9)sU(s 


(3.71) 




4 mb mB mx*{—m\ + m‘j^*)s 1 

mB 


+4 A mB s< 2 mb {%{s) 


2E 


8 Xmb 7)(s)|2(m| — m\*)Ex* + mB s 
T2{s)) - A2(s)(±Cio - Cg){mB - mx*)^ 


+ {mB - mx*){mB + mx*Y{m\ - m%* - s)|4 7)(s) mb Ex*{-mB + mx*) 


+s mB Ai(s)(±C'io — Cg) 


(3.72) 


It is interesting to observe that in the Large Energy Effective Theory (LEET), both helicity am¬ 
plitudes |i7:J^’^(s)p and \H^’^{s)\‘^ have essentially one dependence on the universal form factor 
\ However the helicity amplitude \E[q’^{s)\‘^ is more model-dependent, since it depends on 
the two form factors ^ and \ 

In Eigs. (3.6) and (3.7) we show respectively the helicity amplitudes \H^ _ (-s) P and \H^ _ (s) P 
at leading and at next-to-leading order. We remark that both helicity amplitudes |i7:[^_(s)p are 
completely negligeable comparing to their left-helicity components. Moreover, the impact of the 
NLO corrections on the the amplitudes \H^_{s)\‘^ and \H^_{s)\‘^ increase considerably their 
magnitude up to ~ 100% in the small lepton invariant mass (s < 2 GeV^). The NLO uncertaities 
are dominated mainly by the 5-meson light-cone distribution amplitudes \b,+, the 5-decay con- 
stant /s and the form factor . The corresponding errors were calculated by varing mainly 
these parameters in the indicated range, one at a time, and adding the individual arrors in quadra¬ 


ture. 
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s (GeV2) s (GeV^) 


Figure 3.6: The helicity amplitude (left-handplot) and \Hf(s)\‘^ (right-handplot) at 

NLO (solid center line) and LO (dashed). The band reflects theoretical uncertainties from the 
input parameters [28], 




Figure 3.7: The helicity amplitude (left-handplot) and (right-handplot) at 

NLO (solid center line) and LO (dashed). The band reflects theoretical uncertainties from the 
input parameters [28]. 
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Using the above helieity amplitudes and taking the narrow resonanee limit of K* meson, i.e., 
using the equations 


Tx* — 


9k*Ktt'^K* 


lim 


487r 






= 7i6{p‘^ - , 


(3.73) 


K* 


we ean perform the integration over in Eq. (3.28) and obtain the fourth differential angular 
distribution of i? — K*{-^ Kir) f+^^with respeet to dilepton mass squared s, the azimuthal 
angle 0, the polar angles 9k and 6^+, 


d^T 


3 al^Gl^/Xml\VtbVt*f 


ds d cos 9k d cos 9+ dcj) 128(27r)®m¬ 


(3.74) 


B 


x^^4:Cos‘^9k sin^9+ (^|i7o^(s)p+ |i7^(s)|^j 

+sin^9K (1 + cos^9j^) + \H^[s)\‘^ + \H^[s)\‘^ + |i7:^^(s)|^^ 

-2 sin^9K sin^9+ cos20 Re(^H^{s)H^*{s) + H^{s)H^*{s)^ 

-sin2(t) Im[H^{s)H^*{s) + 

-sin29K sin29+ coscj) Re(^H^{s)H^*{s) + H^{s)H^*{s) + H^{s)H/^*{s) + H^{s)H/^*{s)^ 
-sinct) Im{H^{s)H^*{s) - H^{s)H^*{s) + H^{s)H!^*{s) - 
-2 sin^9K cos9+ (|77^(s)|2 - |77^(s)|2 - \H^{s)\^ + |i7f (s)|2) 

+2 sin9+ sin29K coscj) Re(^H^{s)H^*{s) - H^{s)H^*{s) - Hi^{s)H/^*{s) + H^{s)HI^*{s)^ 
-sincj) Im(H^{s)H^*{s) + H^{s)H^*{s) - H^{s)H^*{s) - H^{s)H^*{s) 


3.3.1 Dalitz distributions 


If we integrate out the angle 9k and 9+ from Eq.(3.74), we get the double 0 angular distribution: 


d^B 
d(l) ds 


a: 


.Gl 


= tb- 


3847r^ 
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nv 


3 \V,>.V„ 
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* 12 


27r 




(3.75) 


-cos20 Re{H^{s)H^*{s) +H^{s)H^*{s)^ + sin2(t) Im{H^{s)H^*{s) + H^{s)H^*{s)^^, 


where tb is the life time of the 5-meson, and the various terms in the expansion above ean be 
speeified , as follows : 


i^o(^)r = ii^o"(^)r+K(^)p, 


(3.76) 
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Figure 3.8: Dalitz distribution 


d^B 


(icos0+ ds 


for B ^ K* i+ r [28] . 


Similarly, we ean get respeetively the 6k and 0+ angular distributions as following: 


d^B 


a 


d cos 9k ds 


= Tb- 


/~i2 

emr-’ F 

3847r® 


^^1 7 hcos^OK \Ho{s)\ 

m'g 4 t 


+sin^dK (^|Fr+(s)|^ + |iF_(s)|^j|, 


(3.77) 


and 


d^B 


d cos 9^ ds 


= Tb- 




3847r' ■ ''™3 
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B 


+ (1 + cos9+y |i7:^(s)|^ + (1 - cos9+)'^ |i7+(s)|^ 
+(1 - cos9+)‘^ |iff (s)l^ + (1 + cos9+)‘^ \H^{s)\'^ 
d^B\H,\^ ^ d^B\H.\^ ^ d^B\H^\2 


(3.78) 

(3.79) 


d cos 0+ ds d cos 0+ ds d cos 0+ ’ 

where the various terms in Eq. (3.79), namely (i^i3|//Q|2/(icos6'+ ds, dfB\ h_\ 2 /dcos,9ds and 
dfB\Hj^\‘ 2 ./d cos 9+ ds can be defined respectively,as follows: 

3 


d cos 6*+ ds 
d%H.\^ 
d cos 6*+ ds 

d cos 6*-i- ds 


= Tb- 


= Tb- 


oLGf 




(3.80) 


^ Tb 


^2 fy2 
-^em^F . 

3847r® 

2 ^2 
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8 ~ cos«+)^ + (1 + cose+y |fl'®(s)|'^|, 

384^5 ^ {(1 + TosOBf |ffi(i>)r + (1 - COBe+f |ff7s)| 
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ds (icos0+ 10 ^ 


Eigure 3.9: Partial Dalitz distribution ^ for B K* £+ £ [28] . 
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Note that the polar angle distribution functions in Eqs. (3.77), (3.78) and (3.80) depend 
only on the modular square terms of the helicity amplitudes, which give the decay width of 
the semileptonic decay (see next section). 

Using our central input parameters given in Tables (3.1) and (A.l) (see Apendix A.l), we 
show in Figs. (3.8), (3.9) and (3.10) at the NLO accuracy the total Dalitz distribution (fB/dcos6+ ds, 
the two angular partial distributions /d cos ds and d‘^B\Ho\‘ 2 / d cos 6*+ ds, respectively. 

From the experimental point of view these Dalitz distribution can serve as a double check of 
whether the branching fraction is different from the SM predictions. 


3.3.2 Dilepton mass spectrum and Forward-backward asymmetry 

Finally, after integrating over the polar angles 0, 9^ and 9k, we derive the total differential 
branching ratio in the scaled dilepton invariant mass for B K* i+i-. 


dB 

ds 


dB\H+\'2 ^ dB\H_\2 ^ dB\HQ\^ 
ds ds ds ’ 


where the various terms in the Eq. (3.81), read as 


(3.81) 


dB\H+\'^ 
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ds 

dB\Ho\2 
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^ 3847r5 
^ 3847r5 

n /2 (^2 

^em'^ F 

^ 3847r5 





\Hjs) 




(3.82) 


The partial lepton invariant mass spectrum p/ ds, dB^n+i^/ds and dB^Uoi^/ds are shown 

in Fig. (3.11) showing in each case the next-to-leading order and leading order results. We remark 
that the partial single distribution dB\H+\' 2 /ds is completely negligeable comparing to the others. 
In fact this is due to the smallness of the helicity amplitude |i7+P, as it is shown in Figs. (3.6) 
and (3.7). 

In Fig. (3.12), we plot the total dilepton invariant mass dB/ds at next-to-leading order and 
leading order. As it is shown in Figs. (3.11-upper plot) and (3.12) the total decay rate is domi¬ 
nated by the contribution of the helicity |i7_ | component. 

We Note that the next-to-leading order correction to the lepton invariant mass spectrum in 
B —is significant in the low dilepton mass region (s < 2 GeV^) but small beyond that 
shown for the anticipated validity of the FEET theory (s < 8 GeV^). Apparently rather large 
uncertainty of our prediction is mainly due to the form factors with their current large uncertainty 
and to a lesser extent respectively due to and the 5-decay constant. 
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ds 10“^ 



ds 10-10 


ds 10 1" 



0.5 

0.4 



0.1 


0 ' ' ' ' 1 ' ' ' ' 2 ' ' ' ' 3 ' ' ' ' 4 ' ' ' ' 5 ' ' ' ' 6 1 

s (GeV2) 


d^\H_ 12 dBr I12 

Eigure 3.11: The dilepton invariant mass distributions —— (upper-plot), — (middle- 

dB\H |2 

plot) and — (lower-plot) for B —atNLO (solid center line) andLO (dashed).The 
band reflects the theoretical uncertainties from input parameters [28,29]. 
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dB 

ds 10-7 



Figure 3.12: The total dilepton invariant mass distribution for B —> at NLO (solid 

center line) and LO (dashed). The band reflects theoretical uncertainties from the input parame¬ 
ters [ 28 ]. 


Besides the differential branehing ratio, B K* i+i- deeay offers other distributions (with 
different eombinations of Wilson eoeffieients) to be measured. An interesting quantity is the 
Forward-Baekward (FB) asymmetry defined in [19,49] 


dvApB 

ds 



dfi 


d^r 

dit ds 


+ 



dfi 


d^r 

dit ds 


(3.83) 


where the variable it eorresponds to 6^, the angle between the momentum of the 5-meson and the 
positively charged lepton £+ in the dilepton CMS frame, through the relation it = —u{s) cos 0+, 
and bounded as 


with 


—fi(s) <u< u(s) , 



(3.84) 


(3.85) 
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2 4 6 

s (GeV2) 


Eigure3.I3: The forward-backward asymmetry in B —> decay at NLO (solid center line) 

and LO (dashed). The band reflects the theoretical uncertainties from the input parameters [28]. 


It is interesting to observe that at the leading order in the LEET approaeh, the EB-asymmetry 
in 5 —deeays depends on one universal form faetor \s), and reads as follows 


dA 


FB 


ds 




2%5 
xCio 


+ ml, + s) + 2 ^ 
s ms 


a 


eff 


7 




(3.86) 


It has been noted in [93] that the loeation of the forward-baekward asymmetry zero Sq is 
nearly independent of partieular form faetor models. An explanation of this faet was given in 
[31], where it has been noted that the form faetor ratios on whieh the asymmetry zero depends 
are predieted free of hadronie uneertainties in the eombined heavy quark and large energy limit. 
Thus the position of the zero sq is given by 


Re(C9"'(l„)) 


’fCr^ 1 1 ^ + 1 

So [ 1 + ml, - s 


(3.87) 


whieh depends on the value of mi, and the ratio of the effeetive eoeffieients Cf /Re (sq) 
Thus, the preeision on the zero-point of the EB-asymmetry in 5 —> is determined 
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essentially by the preeision of the ratio of the effeetive eoeffieients and mi,K We find the in¬ 
sensitivity of So to the deeay form faetors in B ^ a remarkable result, whieh has also 

been diseussed in [93]. However, the LEET-based result in Eq. (3.87) stands theoretieally on 
more rigorous grounds than the arguments based on seanning a number of form faetor models. 
Our result for EBA is shown in Eig. (3.13) to EO and NEO aeeuraey. With the eoeffieients 
given in Table (3.1) and mh = 4.6 GeV, we find the EO loeation of the EB-asymmetry zero is 
So ~ 3.4 GeVl 

In [63] the effeet of the (faetorizable) radiative eorreetions to the form faetor has been studied 
and has been found to shift the position of the asymmetry zero about 5% towards larger values. 
However the effeet of both, faetorizable and non-faetorizable radiative eorreetions modify eon- 
siderably the loeation of the EB-asymmetry zero sq. As it is shown in Eig. (3.13), the numerieal 
effeet of NEO eorreetions amounts to a substantial enhaneement of the EB asymmetry for in¬ 
termediate lepton invariant mass (s = 1.5 — 6 GeV^) and a signifieant shift of the loeation of 
the EB-asymmetry zero to Sq — 4.7 GeV^. The dominant uneertainty (between 5% and up 55%) 
is shared mainly between the 5-meson light-eone distribution amplitudes \b,+, the 5-deeay 
eonstant fs and the form faetor 5 


3.3.3 Transversity Amplitudes for B K*£^£ 

The deeay B —> J/'ipK* is deseribed by three amplitudes** (.4,*; i = 0, ||, ±) in the transversity 
basis, where 4lo(s), 4l||(s) and 41_l(s) have CP eigenvalues -fl, +1 and —1, respeetively [94- 
96]. Here, 4lo(s) corresponds to the longitudinal polarization of the vector meson K* and 4l||(s) 
and 41_l(s) correspond to parallel and transverse polarizations, respectively. The relative phase 
between the parallel (transverse) amplitude and the longitudinal amplitude is given by 0||(±) (s) = 
arg(4l|j(x)Mo(s)). 

The transversity frame is defined as the J/%£ rest frame (see Eig. (3.14)). The K* direction 
defines the negative x axis. The K-k decay plane defines the {x, y) plane, with y oriented such 
that Py{K) > 0. The 2 : axis is the normal to this plane, and the coordinate system is right-handed. 
The transversity angles 6tr and (t)tr are defined as the polar and azimuthal angles of the positively 
charged lepton from the J/ip decay; 6k* is the K* helicity angle defined in the K* rest frame as 
the angle between the K direction and the direction opposite to the J/ip. This basis has been used 
by the CEEO [32], CDE [33], BABAR [34], and the BEEEE [35] collaborations to project out the 
amplitudes in the decay B J/ipK* with well-defined CP eigenvalues in their measurements 
of the quantity sin 2{3, where /d is an inner angle of the unitarity triangle. 

We also adopt this basis and analyze the various amplitudes from the non-resonant (equiv¬ 
alently short-distance) decay B K*Bi~. In this basis, both the resonant B K*J/ip —> 

I'the corresponding quantity in the inclusive decays B , for which the zero-point is given by the 

solution of the equation Re^C 9 ®®(so)^ = 

**they should not be confused with the form factors Ao(s), Ai (s) etc. 
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Eigure 3.14: Definitions of the transversity angles 9tr, fitr, cind 9k*- The angles 9tr and fitr are 
determined in the J/fi rest frame. The angle 9 k* is determined in the K* rest frame. 


(already measured) and the non-resonant {B amplitudes turn out to be very 

similar, as we show here. 

The angular distribution is given in terms of the linear polarization basis (.4,±i (s) = (s) ± 

A±{s))/^/2) and ^o(s) by 


__ 

ds d cos 9tr d cos 9 k* dfitr 


ffiw) ■ +/2H ■ iA(^)r+ fsH ■ i^^(^)r 

riffiw) ■ Im(^|(s)^_L(s)) 4- fi^w) ■ Re(^o('5)‘^l|('S)) 
rjfeiw) ■ Im(^o(s)^x(s)) , 


where r] = +1(—1) for B^ and (B^ and B ), and the coefficients /j=i,...,6, which depend on 
the transversity angles w = {9K*-,9tr-, fitr), are given by: 


fl{w) = 9/{327 i) ■ 2 COS^ 9k* {1 — silfi 9tr cos'^ fitr), 

f 2 {w) = 9/{327r) ■ sin^ 9 k*— sin^ 9trSin^ fitr), 

fsiw) = 9/{327 i) ■ sm^ 9 k* sin^ 9tr, 

ffiw) = 9/{32Tr) ■ sin^ 9 k* sm29tr sin fitr, 

U{w) = - 9/(327r) ■ I/V 2 ■ sin 29 k * sin^ 9tr sin 2(j)tr, 
hiw) = 9/(327r) ■ I/V 2 ■ sin 29 k * sin 29tr cos fitr ■ 
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Group 

O 

to 


I^IlP 

0± 

011 

CEEO [32] 

0.52 ±0.08 

0.16 ±0.09 

0.32 ±0.12 

-3.03 ±0.46 

-3.00 ±0.37 

CDF [33] 

0.59 ±0.06 

n 1 q-i-0.13 

0.28 ±0.12 

-2.58 ±0.54 

-2.20 ±0.47 

BaBar [34] 

0.60 ±0.04 

0.16 ±0.03 

0.24 ±0.04 

-2.97 ±0.17 

-2.50 ±0.22 

Belle [35] 

0.60 ±0.05 

0.19 ±0.06 

0.21 ±0.08 

-3.15 ±0.21 

-2.86 ±0.25 

AS [28] 

0.51 

0.21 

0.28 

-3.25 

-3.04 


Table 3.3: Current measurements of the decay amplitudes in the transversity basis for the decay 
B — J/fK*. The corresponding amplitudes for the non-resonant decay B worked 

out in this paper in the LO approximation at given in the last row [28], 


In terms of the helicity amplitudes introduced earlier, the amplitudes in the linear polar¬ 

ization basis, ^o,±,||> can be calculated from the relation: 


A ( s ) = k{h^{s) + , 

^±i(s) = K{Hi{s) + Hi{s)) , 


with 

3847r^ rriQ 


|viw, 


* 12 
ts \ 


Experimental results are conventionally expressed in terms of the spin amplitudes 4,o,±,|| 
normalized to unity, with + 1*4111^ = 1. We show the polarization fractions, 

Eo/r = |4.o(s)p, Eii/r = |4.||(s)p and r^/E = |4 ._l(s)P in the leading and next-to-leading 
order for the decay B in Fig. (3.15), respectively. Since the interference terms in the 

angular distribution are limited to Re(^ll^o), Im(^x^g) and Im(^x^j|), there exists a phase 
ambiguity: 


-011 , (3.88) 

±7r - 0x , (3.89) 

dzTT - (0_L - 0||) . (3.90) 

To avoid this, we have plotted in Fig. (3.16) the functions cos0||_j_(s) and sin^y j^(s), showing 
their behaviour at the leading and next-to-leading order. The dashed lines in these figures corre¬ 
spond to using the FO amplitudes, calculated in the FEET approach. In this order, the bulk of the 
parametric uncertainty resulting from the form factors cancels. Although, strictly speaking, the 
domain of validity of the FEET-based distributions is limited by the requirement of large energy 
of the K* (which we have translated into approximately s < 8 GeV^), we show this distribution 
for the entire s-region allowed kinematically in 5 —> The shaded curves correspond to 

using the NFO contributions in the FEET approach. 


0± 


0± - 
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Eigure 3.15: The helicity amplitudes |^o('S)P (upper-plot), |^||(s)p (middle-plot) and |^±(s)P 
(lower-plot) in B at NLO (solid center line) and LO (dashed). The band for NLO 

reflects theoretical uncertainties from input parameters [28]. 
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Figure 3.16: The functions cos0||^x(s) and sin0|| x(s) at NLO (solid center line) and LO 
(dashed). The band reflects all theoretical uncertainties from parameters with most of the uncer¬ 
tainty due to the form factors \D).The vertical line at s = 8 GeV^ represents the domain of 

validity of the LEET approach in our case [28]. 


We eompare the resulting amplitudes and 0±(s) at the value 

s = with the eorresponding results from the four experiments in Table (3.3). In eompar- 
ing these results for the phases, we had to make a ehoiee between the two phase conventions 
shown in Eq. (3.90) and the phases shown in the last row of this table correspond to adopting 
the lower signs in these equations. We note that the short-distance amplitudes from the decay 
B —> are similar to their resonant counterparts measured in the decay B —> J/flK*. 

We also note that a helicity analysis of the decay B —> J/flK* has been performed in the QCD 
factorization approach by Cheng et al. [97]. 

The structures in the phases shown in Fig. (3.16) deserve a closer look. We note that at the 
leading order, the phases (p±{s) and 0|| (s) are given by the following expressions: 


0±(s) 

0||(s) 


Arg 

Arg 


—is + 2mi,mB | (g)] - Arg[^o(s)] , (3.91) 

'mbmB ys t J J 

^(5) f^ eff , o „„ „„ ^ eff 


rriby/s 


s Cfl + 2 mf, C 7 ' 


-2mbmB (+ 


2mbmB 


,m 




m 


B 


Arg[ylo(s)], 


(3.92) 
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Eigure 3.17: The phase 0o('S) at NLO (solid line) and LO (dashed) [28]. 


where we can neglect the term proportional to (my/m?^) in the latter equation. The phase 
0o(s) = Arg[^o(’S)] is constant in the entire phase space, as shown in Eig. (3.17). The func¬ 
tions in the square brackets in Eqs. (3.91) and (3.92) are purely imaginary. However, due to 
the fact that in the SM the coefficients Cg® and have opposite signs, these phases become 
zero at a definite value of s, beyond which they change sign, yielding a step-function behaviour, 
shown by the dotted curves in the functions cos0||(s) and cos(j)±{s) in Eig. (3.16), respectively. 
The position of the zero of the two functions, denoted, respectively, by Sg and sjj, are given by 
solving the following equations: 


Arg 


i \/X 


nibmB V-Sq 
-iEK> 


Arg 


Sg C9'^®(so)-f ^(sg) 

JC9®®(s|j) -f ^(s|j) 


mb V55 


00 (Sq) ) 

(3.93) 

00 (4)- 

(3.94) 


Eor the assumed values of the Wislon coefficients and other parameters, the zeroes of the two 
functions, namely Sg and Sg , occur at around s ~ 3 GeV^, in the lowest order, as can be seen in 
Eigs. (3.16), respectively. 

The LO contributions in sin0||(s) and sin0x(s) are constant, with a value around 0, with 
a small structure around s ~ 3GeV^, reflecting the sign flip of the imaginary part in ^||(<s) 
(•A_l(s)). At the NLO, the phases are influenced by the explicit 0 { as ) contributions from the 
factorizable and non-factorizable QCD corrections (see section 3), which also bring in parametric 
uncertainties with them. The most important effect is that the zeroes of the phases as shown for 
cos0_l(s) and cos0||(s) are shifted to the right, and the step-function type bahaviour of these 
phases in the LO gets a non-trivial shape. Note that in both figures a shoulder around s ~ 8 GeV^ 
reflects charm production whose threshold lies at s = 4m^. 
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3.4 Decay Distributions in B ^ p 


After a complete analysis of the B we turn now to the semileptonic B ^ p iui 

one. In this section, we present general spectra analysis in exclusive B —> piu^ decay in terms 
of the corresponding helicity amplitudes. Further, we calculate the different dilepton invariant 
mass distributions and the corresponding Dalitz distributions. We include the 0(Q;s)-corrections, 
1/E power corrections by means of the large energy expansion technique (LEFT) and using the 
light-cone QCD sum rules approach. 

First, let us describe the apropriate matrix elements forB^p iui. Since this decay is purely 
di iy — A) transition, one could get the corresponding matrix element from the B 
ones, using the following replacements : 


Co = 


-Cio - 2 , 


G p C(p 


V2 


71 


G 




eff 


= 0 , 


4Gf 


■K 


ub 


(3.95) 

(3.96) 

(3.97) 


This amounts to keeping only the charged current (V — A) contribution m B ^ K*Gi~ decay. 
Thus, the corresponding amplitude for the semileptonic b ^ u decay, can be factorized into 
a leptonic and a partonic part as, 

M{h uiue) = -i^Vub^[u'y^{l - 75 ) b] [1Y{1 - 75 )^'<;]}- (3-98) 

From the semileptonic amplitude given in Eq. (3.98), we notice that the exclusive b ^ u 
decays is a good candidate for a clean determination of the modulus of Vub, one of the smallest 
and least well known CKM matrix elements. Experimentally, the main difficulty of the obser¬ 
vation of 6 —> signal events is the large background from b cfz/^ events'!”^. For that, 

different experimental distribution analysis are in order to overcome this prblem. In this spirit, 
we propose many angular distributions studies of B ^ p fz/^, where the vector meson decays to 
two pseudoscalars^^ p tt+tt” . 

Four independent kinematic variables completely describe the semileptonic decay B ^ py 
7r~^7r~) y the four variables most commonly used are the invariant dilepton mass distributions 
and three polar angles. Thus, the differential decay rate for B p(— 7r^7r~)iui is expressed in 
terms of the helicity amplitudes H± o by [98]: 


—77 = —77—77GplVubl — p-Bip tt^tt ) 

ds d cosOp d cos 6d(j) 8(47r)"^ m\ 

ttBecause \Vub/Vcb\ ~ 0.1, the branching fractions of the exclusive b uGi decays (~ 10“^) are small 

compared to those of the charmed semileptonic decays, which are of the order of some percent. 

Mthis is due to the fact that the non-leptonic p tt+tt” decay is by far the dominant branching ratio [36], with 

B{p 7r+7r“) ~ 100%. 
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X 1(1“ cosO+Ysin^Op 

+ (1 + cos6+Ysin^9p |iJ_(s)|^ 

+4 siv?9+ cos^9p |Ero(s)|^ 

—4 sin9-^. (1 — cos9+) sin9p cos9p coscj) H+{s)Ho{s) 

+4 sin9+ (1 + cos9^) sin9p cos9p cos(p H_{s)Ho{s) 

—2 sin‘^9^ sin^9p cos2<p H^{s)H_{s)^. (3.99) 


The function A can be found in subsection 3.1.2. The angles here 9^, 9p and 0 are defined re¬ 
spectively as : the direction between the charged lepton and the recoiling vector meson measured 
in the W rest frame, the polar angle between 7 r+ (or tt”) and the direction of the vector meson 
in the parent meson’s rest frame, and the azimuthal angle between the planes of the two decays 

B —> p(— 7r~^7r~)iui. 

The helicity amplitudes can in turn be related to the two axial-vector form factors, Ai(s) and 
^ 2 ( 5 ), and the vector form factor, V(s), which appear in the hadronic current [98]: 


H_{s) = {niB + mp) Ai{s) + 2 


Vx 


rriB + m. 


V{s), 


V\ 


H+{s) = {tub + nip) Ai{s) - 2 -^-U(s), 

niB + nip 

1 


HAs) = 


2mp^/s L 


(m^ -m - s){mB + nip) Ai(s) - 4 


A 


A2{s) 


(3.100) 

(3.101) 
,(3.102) 


niB + nip 

where nip stands for the p-meson mass. Using Eqs. (2.62), (2.63) and (2.64) in Eqs. (3.100)- 
(3.102), we obtain at the LO accuracy in the large Energy Limit, the following helicity ampli¬ 
tudes: 


H(s) = 2 


H^{s) 

Hois) 


= 2 


En + 


Ep- 


niB 

tub 


A\>), 


niB nip ^/s L 


tub Ep (m| -ml- s) 


-2 A 




2 A .(p) 


mBEp^/s 


C(^), 


(3.103) 

(3.104) 


(3.105) 


where Ep represents the p-meson energy, defined in Eq. (2.50) (see subsection 2.4.1), and 
denote the two universal form factors for the B piu£ transition in the LEET theory. 

As this framework does not predict the corresponding decay form factors, they have to be 
supplied from outside. Eor that we have suggested the following: 

• Having at hand the appropriate B LEET form factors, namely ^ and \ 

one can relate them easily to the B piug ones in the SU(3)-symmetry limit. Eollowing 
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Figure 3.18: LEET form factors ||('5) in B ^ p Ivf,. The central values are represented by the 

dashed curves and the bands reflect the uncertainties on the form factors [28], 


this statement, then the semileptonic B piug LEET form faetors ean be defined as 
'( 0 ) = «Hii (0). Unhappily, one has to eonsider the SU(3)-breaking effects in the 
corresponding form factors, which have been evaluated within the QCD sum-rules [99]. 


Thus, we take the SU(3)-symmetry breaking factor, as 


(3.106) 


Taking this and uO) from Table (A.l) into account, we obtain the corresponding one 
for the B p-transition: 

= 0.22 ±0.04. (3.107) 


• To extrapolate the B —> piu^^ form factors at s 7 ^ 0, we use the same extrapolation function 
as for B —> form factors (see Eq. (3.69) in subsection 3.2.3): 



(3.108) 


Using Eq. (3.108), we have plotted in Eig. (3.18) the corresponding B ^ p EEET from factors, 
namely and ^[|^^(s). To check the consistency of the corresponding form factors, we have 

compared the behavior of .^^|| (s), with the one used by [30] and surprisingly it turns out that the 
agreement is reasonable. 
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Eigure 3.19: Dalitz distribution 


d^B ^ „ 

dcos6*+ 


piiye [28]. 


We notice that the apparently rather large uncertainty of our prediction, typically ~ ±25%, is 
mainly due the the B —> K* form factor ■Ci || (0) form factors with their current large uncertainty 
and to a lesser extent due to the SU(3)-symmetry breaking factor Csu{ 3 )- It may be hoped that in 
the longer term future the form factors could be known with much greater confidence. 


3.4.1 Dalitz distributions 

Integrating out the angles {9p, 9+) and (0, 9+) we obtain respectively the double 0 and 9p angular 
distributions: 


and 


d^B G2 

d^s - 

\Ho{s)\^ + \H_{s)\^ + \H+{s)\^ - cos2(j) H_{s) H+{s)], (3.109) 


d^B 


d cos 9p ds 




2 cos^9p \Ho{s)\‘^ + sin^9p ^|if+(s)|^ ± |if-(s)|' 


(3.110) 


Similarly, We give here the 0+ double angular distribution as following: 
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Figure 3.20: Dalitz distribution 


dcos^p 


pii^i [28]. 


d^B 

d cos 6*+ ds 


tb 


^2 
Or p 


s/X 


256mg7r^ 


\Vub\‘^{B{p ^ Tl^Tl )) 


|2 sin^Qj^ |ifo(s)P + (1 — cosOj^Y I-^+('5)P + (1 + cosOYf |if_(s)p| 

d%n,\^ ^ d%u,\^ ^ d%H.p ^ 

d cos 0+ ds d cos 6*+ ds d cos 0+ ds ’ 


where the corresponding partial double angular distributions read as: 


d cos 0+ ds 

d%H.\^ 
d cos 0+ ds 
d%H,\^ 
d cos 0+ ds 


^ - cosO+Y |i^+(s)|'} 

^ 7r+7r-))|(l + cosO+Y l^-(s)P} 

Tb IOfer iB{p 7r+7r“))|2 sin^e+ |ifo(s)p|. (3.112) 


Implementing the 0(as)-improvements in the various helicity amplitudes above, we have 
shown respectively in Figs. (3.20) and (3.19) the explicit behavior of the 9p and 9+ angular distri¬ 
butions. Wheras in Figs. (3.21) and (3.22) we have presented the 6*+ partial angular distributions: 

((i^i3|j:^_|2/(icos6*+ ds) and ((i^i3|j:/Q|2/(icos6'+ ds). 
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(fiB, 


H_ 


Eigure 3.21: Partial Dalitz distribution ^ for B piu^ [28]. 



dfB,^ |2 

Eigure 3.22: Partial Dalitz distribution ds ^^^ ^ 
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To project out experimentally the various helicity components namely ITTqP and one 

can use the 9p and 6*+ Dalitz distributions respectively. Extracting the ones, which is 

the more easiest one**, requires a precise measurement at = 0. This means that the charged 
lepton should be back to back to the recoiling p-meson. However measuring iTToP imply a higher 
efficiency at = 0 or tt, by means that the vr-meson should be measured in the same axe as the 
recoiling p-meson. 


3.4.2 Dilepton mass spectrum 

Finally, integrating out the polar angle 6*+, 9p and (p from Eq. (3.99), we obtain the total branching 
decay rate for the B piu^ transition: 


ds 


^ tt+tt ))||iJo(s)P + + |i/-(s)p} 

dB\HQ\2 ^ dB\H+\-2 ^ dB\H_\‘2 
ds ds ds 


(3.113) 


The partial distribution amplitudes defined in Eq. (3.113), read as follow: 


dB 


|rr_|2 


ds 

dBiH+i'i 

ds 

dB\Ho\^ 


ds 


= Tb 


= Tb 


= Tb 


Gls^/X . 

96m^7r^ 

GlsVX 
96m^7r^ 
GlsVX^ 
96ml 


B 




\\B{p- 

7r’''7r“))<| 


\\B{p- 

7r’'"7r“))<| 


\\B{p- 

7r’'"7r“))<| 

[|rro(s)f} 


(3.114) 


In Figs. (3.23) and (3.24), we have plotted respectively the various partial dilepton invariant 
mass distributions and the total one. 

The observed helicity \H^ \ componenE^ in Fig. (3.23-middle plot), is completeley neglige- 
able comparing to the two others helicities. This result is just a direct consequence of the {V — A) 
coupling. The semileptonic process A b ^ ul V produce a u quark which is predominately he¬ 
licity A = —1/2. In a i? —process, the helicity of the X meson is then determined by 
whether the u quark combines with a spectator quark that has helicity A = -1-1/2 or A = —1/2. 
If X is a spin-zero meson, only A = -1-1/2 spectator quark contributes. However, if X has spin 
1, both helicities of the spectator quark contribute, leading to X helicities of A = 0 and A = — 1, 
but not A = -fl. 

** this is due to the fact that |iJ_ | « see Eq. (3.100). 

'1'° Note that from Eq. (3.114)—« |iJ_i_(s)p. 

t * the same argument holds for b ^ c process. 
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d^\H^\2 
ds |U^b|2lO-4 




dB 




dBijj ,2 

plot) and — (lower-plot) for B 


ds 


(middle- 


piof at NLO (solid center line) and LO (dashed).The 


band reflects the theoretical uncertainties from input parameters [28,29]. 
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dB 

ds \V^b\‘^ 



Figure 3.24: The total dilepton invariant mass distribution for B —> piui at NLO (solid center 
line) and LO (dashed). The band reflects theoretical uncertainties from input parameters [28]. 


Consequently, the B p i~U£ transition should be completely dominated by the two helicity 
components \Hq\ and \H_\, in a good agreement with what we have observed in Fig.(3.23 upper- 
plot and lower-one). Contrary to the B decay rate, the B —> piui one is totally 

dominated by the helicity \Hq\ component, as it is shown Figs. (3.23-lower plot) and (3.24). 

Note that the impact of the NLO correction on the total branching ratio and the partial ones 
are less significant as for the B decay, and these is simply due to the absent of the 

penguin form factor corrections in the B piue decay. However the large systematic error of 
our prediction comes from the uncertainty in the form factors with their current large uncertainty 
and to the SU(3)-breaking effects. There is no doubt that a precise measurement of the long¬ 
distance effect, will reduce considerably our uncertainty on B ^ piui decay. 


3.5 Phenomenological Discussion on Rb 

In the Standard Model, the charged current weak interactions of three generations of quarks 
are governed by a Lagrangian which contains a transformation from the mass eigenbasis to the 
flavour (generation) eigenbasis [36,100,101]. This flavor-mixing is expressed as a 3x3 complex 
matrix V^km known as the Cabibbo-Kobayashi-Maskawa (CKM) matrix [5] (see Eq. (2.4)). 
The Unitarity of this matrix reduces the number of independent parameters to nine, which 
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can be chosen as three real mixing angles and six imaginary phases. Live of the phases are 
removable. The four remaining parameters are fundamental constants of nature, to be determined 
by experiment since the SM itself gives no guidance as to their values. Therefeore, an important 
target of particle physics is the determination of the CKM matrix [5]. 

Eig. (3.25) illustrates the hierarchy of the strengths of the quark transitions mediated through 
charged-current interactions^^, and Table (3.4) their present experimental source of informa¬ 
tion [102]. 



- 0 ( 1 ) 

0(10 S 

--- O(io") 

0(10 ^ 


q=- 1/3 Q=-r2/3 


Eigure 3.25: Hierarchy of the quark transitions mediated through charged currents [103]. 


It it known that the only CP violation source in the SM is supposed to arise from a single 
phase in the CKM matrix. This is a very remarkable property of the Kobayashi-Maskawa picture 
of CP violation: quark mixing and CP violation are closely related to each other. This property 
is often used to determine the angles of the unitarity triangle without the study of CP violating 
quantities. 



Eigure 3.26: Unitarity triangle. 


The measurement of Rh = \VublVcb\ (= \Vub/VtbVff\) constraints the length Ru of the unitar¬ 
ity triangle (see Eig .(3.26)) through the relation [100]: 


R 


U 


v:bVud 

VfbVc, 


(3.115) 


Transitions within the same generation are governed by CKM elements of 0(1), those between the first and the 
second generation are suppressed by CKM factors of 0(10“^), those between the second and the third generation 
are suppressed by 0(10“^), and the transitions between the first and the third generation are even suppressed by 
CKM factors of 0(10-3). 
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Relative amplitude 

Transition 

Source of information 

~ 1 

u ^ d 

Nuelear /9-deoay 

~ 1 

c ^ s 

Charmed partiele deeays 

~ 0.22 

u ^ s 

Strange partiele deeays 

~ 0.22 

c ^ d 

Neutrino eharm prod. 

~ 0.04 

c ^ b 

b deeays 

~ 0.003 - 0.004 

u ^ b 

Charmless b deeays 

~ 1 

t ^ b 

Dominanee of t ^ Wb 

~ 0.04 

t ^ s 

Only indireet evidenee 

~ 0.01 

t ^ d 

Only indireet evidenee 


Table 3.4: Relative strengths of charge-changing weak transitions [102]. 


While the two elements Vud and Vcd are known with high aeouraey [36,100], the two left in Ru, 
namely \ Vub\ and | Kb |, are under extensive diseussion at present, espeeially |Kfe|- Their values are 
measured mainly in semileptonie 5-deeays using two independent methods. An endpoint anal¬ 
ysis in inelusive semileptonie 5-deeays yields a direet determination of |14fe/V[;b| [104], while 
measurements of branehing fraetions of exelusive final states sueh as B ^ (tt, p)fz/£ measure 
iKiftI [105]. The model-dependenee in either method is quite substantial. Sinee we have ana¬ 
lyzed the exelusive deeays B — K*£^£~ and B — p£o£, we will propose a model-independent 
analysis of this ratio. Before doing that, let’s have some inside about their experimental status. 


3.5.1 Rb Phenomenology 

Exelusive semileptonie b u£vi deeays are an aetive area of experimental and theoretieal 
study [106-127]. These rare proeesses ean be used to extraet the magnitude of Vub, one of 
the smallest and least well known elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark¬ 
mixing matrix [5]. Beeause iKb/Kfel ~ 0.08, the branehing fraetions for exelusive b u£iy£ 
proeesses are small, of order 10“^, and they have only reeently beeome experimentally aeeessi- 
ble. 

Extraeting |14fe| from a measured deeay rate requires signifieant theoretieal input beeause 
the matrix elements for sueh proeesses involve eomplex strong-interaetion dynamies. Although 
the underlying b u£iye deeay is a relatively simple weak proeess, it is diffieult to ealeulate 
the strong-interaetion effeets involved in the transition from the heavy B meson to the light 
daughter meson. Beeause of these theoretieal uneertainties, even a perfeetly measured B p£vi 
branehing fraetion would not at present lead to a preeise value of \ Vub\- 

The dynamies in B ^ p£vi deeay are in eontrast with b c£ui deeays, sueh as B D*£ui, 
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where a heavy quark is present both in initial and final states. In this case, techniques based 
on HQET can be used to calculate the decay amplitude with good precision, particularly for the 
kinematic configuration in which the charm hadron has zero recoil velocity. The zero-recoil point 
in B ^ piu( cannot be treated with similar techniques, however, because the daughter u quark is 
not heavy compared to the scale of hadronic energy transfers. Nevertheless, substantial progress 
has been made using a variety of theoretical methods, including quark models [111-117], lattice 
QCD [118-120], QCD sum rules [121-122], and models relating form factors measured in 77 —> 
K*iiy£ decay [128] to those in 77 — piu^ decay. 

Experimentally, the main difficulty in observing signals from b uiu^ processes is the very 
large background due to 6 — c£ui. Because a significant fraction of 77 — piu^ events have 
lepton energy beyond the endpoint for b ciu^ decay, lepton-energy requirements provide a 
powerful tool for background suppression. However, extrapolation of the decay rate measured 
in this portion of phase space to the full rate again requires the use of theoretical models, and 
it introduces model dependence beyond that associated with simply extracting the value of | Vub \ 
from the branching fraction. 

The BABAR collaboration has recently presented a preliminary measurement of the CKM 
matrix elements \Vub\ with the charmless exclusive semileptonic 77° —p“e+z/e decay. Their 
result is [129]: 

13(5° ^ p-e+z/e) = (3.39 ± 0.44 ± 0.52 ± 0.60) X 10“^ 

\Vub\ = (3.69 ± 0.23 ± 0.27+°;^°) X 10“°, (3.116) 

where the quoted errors are statistical, systematic, and theoretical respectively. To extract \ Vub\, 
they have used different form-factor calculation. 

In order to reduce the large theoretical errors on the form-factor, it is more convenient to 
study the distribution of s is reflected in the p momentum spectrum. Eventually, studies of the 
s distribution, as well as of the angular distributions of the decay products, should reduce the 
model dependence on \ Vub\ by constraining theoretical models for the decay form factors. In the 
next subsection, we propose a model-independent analysis of the ratio Rb using the s Helicity 
distribution. 


3.5.2 Model-independent analysis of Rb 


To reduce the non-perturbative uncertainty in the extraction of |14fe|, we propose to study the 
ratios of the differential decay rates in 5 —pfz/^ and B K*Ri~ involving definite helicity 
states. These s-dependent ratios Ri{s), {i = 0, —1, -fl) are defined as follows: 

/ ds 


Riis) = 




/ ds 


Erom Eqs. (3.81) and (3.113), one obtains straightforwardly: 


Ri(s) = 


1 77 ) 


(ir*), 




47rs Rl \h['^\s)\^ ’ 


(3.117) 


(3.118) 
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Figure 3.27: The Ratio R-{s) with three indicated values of the CKM ratio Rb = |i4b|/| |. 

The bands reflect the theoretical uncertainty from Csu{ 3 ) = 1-3 ± 0.06 and ^(0) = 0.28 ± 
0.04 [28]. 


where Rb = and the term ^(s) should be understood as the heliety 

amplitudes for the deeay B K*{-^ iF + vr) {B p(—> vr+vr”) [vfj, defined in 

seetion 3.3 (3.4). 

The ratio R- (s) suggests itself as the most interesting one, as the form faetor dependence 
essentially cancels. It is easy to see this cancellation especially in the LEFT approach as the 
function R_ (s) scales in term of form factor as: 

1 ) 2 

^-(®) 0^(^TTTv) (3.119) 

^ C± [S) ^ 

where the function [F{s,mK*,rnp ,...) denotes the kinematical contribution extracted from the 
ratio \H-^ ^(s)p/|i7f^^(s)p. Since the heliety amplitudes \s) and H^\s) depend just on 
one universal form factor, respectively ^ and \ one can easily factorize apriori the ratio 
R-{s) as a kinematic contribution (expressed in R{s, mK*,mp ,...)) and a dynamical contribu¬ 
tion (encoded in the ratio of V'C±^)- 

There is no doubt that in the SU(3) symmetry limit, the function R-{s) will be defined as 
a kinematical function and thus no uncertainty from the non-perturbative regime. Unhappily, 
the reality is far from being that, and one has to incorporate the SU(3)-breaking effects. Then 
the only source of uncertainty coming from the long distance contribution, will be translated 
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Eigure 3.28: The Ratio Ro{s) with three indicated values of the CKM ratio Rb = |I4b|/| |. 

The bands reflect the theoretical uncertainty from Csu{ 3 ) = 1-3 ± 0.06 and ^(0) = 0.28 ± 
0.04 [28]. 


in the ratio Csu{ 3 ) = which turns out to be Csu{ 3 ) = 1-3 ± 0.06, already defined in 

section 3.4. 

On the other hand, the ratio Ro{s) could bring a certain hint on the the structure of Rb. 
However its dependence on the form factors is more involved and cann’t be fudged away, as one 
can see: 


1 I'^l^ Hs) \s) \s,mK*,...)\^ 


(3.120) 


where \s,mK*, ■■■) (s, •••)) ^ certain dynamical function^' obtained from 

|/7(f*(s)p (|i7o('S)P)- Thus, we see that the ratio Rois) is less attractive than the i?_(s) ones for 
the extraction of the CKM ratio Rb, since its uncertainty is more poluated with the long distance 
contribution. 

In Eigs. (3.27) and (3.28) we plot R-{s) and Ro{s) respectively, for three representative 
values of the CKM ratio Rb = 0.08, 0.094, and 0.11. However, as we noticed earlier, the ratio 
R- (s) may be statistically limited due to the dominance of the decay B ^ piu^ by the Helicity- 
0 component. Eor the LEET form factors used here, the compounded theoretical uncertainty is 

To extract the functions ^ , one has just to write \Hf (s)p defined in Eqs. (3.76) 

(3.102), as (3.120) and then identify them. 
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shown by the shaded regions. This figure suggests that high statistics experiments may be able 
to determine the CKM-ratio from measuring Ro{s) at a competitive level compared to the other 
methods en vogue in experimental studies. 


3.6 Summary and Outlook 

In this Chapter we have investigated an 0(Q;s)-improved analysis of the various helicity compo¬ 
nents in the decays B and B piui, carried out in the context of the Large-Energy- 

Effective-Theory. Using that and borrowing the corresponding form factors from QCD sum rule, 
we have investigated the corresponding distributions, decay rates and Eorward-bckward asym¬ 
metry IL 

In the first part, we have concentrated mainly on the study of double and single angular 
distributions, and the EB asymmetry. Our findings can be summarized as follows [28]: 

• We have calculated the helicity components implementing the (^(q;^) corrections and shown 
that the -fl-helicity component is completely negligeable compared to the two other com¬ 
ponents, namely Ho and H_. 

• The total dilepton invariant mass distribution ^ is dominated by the partial single distri- 

bution —^—. The nex-to-leading order correction to the total dilepton invariant mass 
distribution is significant in the low dilepton mass region (s < 2(GeV2)). 

• We have shown that the 0(0;*) effects on the forward-backward asymmetry shifts the pre¬ 
dicted location of its zero by ~ 1 GeV^, confirming essentially the earlier work of Beneke, 
Feldmann and Seidel [25]. 

• We have carried out the B decay analysis in the so-called transversity ba¬ 

sis. We have compared the EEET-based amplitudes in this basis with the data currently 
available on B ^ K* J/ip and find that the short-distance based transversity 
amplitudes are very similar to their long-distance counterparts. 

In the same spirit we have studied the B piui decay in the second part of this chapter, using 
the helicity analysis in the large energy effective theory . After presenting various double and 
single angular distributions, we summarize [28]: 

• Considering the SU(3)-breaking effects, we have related the B piu^ EEET-form fac¬ 
tors, namely and to the corresponding form factors in B ^ . Our 

numerical estimates on ^ (s) and (s) are in agreement with the ones worked out for 
the full QCD form factors in the QCD sum-rule approach in [30]. 

Ifihe FB asymmetry is investigated just for the B K*H£~ decay. 
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• Implementing the 0{as) corrections to the B piui helicity components, we have shown 
that the +l-helicity component is completely negligeable compared to the two other com¬ 
ponents, namely Ho and H_. 

(its 

• The total dilepton invariant mass distribution ^ is dominated by the partial single dis- 

dB\H |2 

tribution —The nex-to-leading order correction to the total dilepton invariant mass 
distribution is completely negligeable. 

Einally, combining the analysis of the decay modes B and B we have 

shown that the ratios of differential decay rates involving definite helicity states, R-{s) and 
Ro{s), can be used for extracting the CKM matrix elements iKfel/lUtsI in a model-independent 
way. 




Chapter 4 


Exclusive B Decay in SUSY 

This chapter is devoted to the semileptonic rare B —> decay, by contrasting its antic¬ 

ipated phenomenological profile in some variants of supersymmetric models. We discuss the 
constraints on the Wilson coefficients Cj, Cg, Cg and Cio, that the current data on rare B decays 
implies in the context of minimal flavour violating model and in more general scenarios admitting 
additional flavour changing mechanisms. As probes of new physics effects in B ^ 
we propose to study the ratios Ro{s) and i?_ (s) (introduced in the previous chapter) using some 
generic SUSY effects. 


4.1 Introduction 

Although the Standard Model (SM) of the elementary particle physics is successful in explaining 
almost all experimental results, it is possible that physics beyond the SM exists just above the 
presently available energy scale. Since new physics may affect various processes at low energy 
such as the flavor changing neutral current (FCNC) processes of iF-mesons and 5-mesons, new 
physics searches in these processes are as important as direct particle searches at collider experi¬ 
ments. A prime example is the 6 —s 7, process. Experimentally the current world average based 
on the improved measurements by the BABAR [11], CLEO [12], AEEPH [13] and BEEEE [14] 
collaborations, B{B Xg'j) = ( 3 . 43 )'^ q 37 ) x 10“"^. It is known that this process puts very strong 
constraints on various new physics beyond the SM, for example two Higgs doublet model and 
supersymmetric (SUSY) extension of the SM. 

Along with the b ^ s j process, another important rare b decay process is the b s i+r 
decay. In particular, after the first measurements of the semileptonic rare 5-decays reported in the 
inclusive 5 —mode, with 5(5 — Xsi~^i~) = (6.1 ± 1.4+};^) x 10“®, 
by the BEEEE collaboration [14] as well as the exclusive 5 —= e^,/i^) one, 
typically 5(5 —= (l.GS^'^oJg ± 0.28) x 10“®, by the BABAR collaboration [18]. 

With increased statistical power of experiments at the 5-factories in the next several years, 
the decays discussed above and related rare 5 decays will be measured very precisely. On the 
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theoretical side, impressive progress in the theoretical precision has been achieved concerning 
the exclusive as well as the inclusive semileptonic (and radiative) rare -B-decays, with the com¬ 
pletion of NNLO (NLO) QCD calculations [21-25]. Although the theoretical uncertainties can 
be addressed only with a complete NLO, as 6 —sy is calculated so far only in NLO , the 
SM value for the branching ratio is in agreement with the experimental measurement within the 
1 — 2a level. 

Since the abovementioned ECNC rare 5-processes are forbidden in the Bom approximation 
any significant deviation from the SM would imply strongly the existence of new physics, such 
as Supersymmetry (SUSY). The reason is simply that, while in the SM the b ^ s transition 
are dominated by one loop contributions with the exchange of a virtual W and the top quark, in 
SUSY [130,131] several competing sources of ECNC are present. To begin with, in SUSY mod¬ 
els the Higgs sector is richer than in the SM, since at least two Higgs doublets must be present. 
Consequently, there exists at least one physical charged scalar which can be exchanged in the 
one-loop contribution to 6 —s, together with an up quark. The second obvious source of ECNC 
comes from the supersymmetrization of the W and the charged Higgs contributions, where the 
up quark is replaced by an up squark and W. 

In this chapter we present a SUSY analysis of the semileptonic rare B —> decay, 

in the so-called minimal supersymmetric standard model (MSSM). After a brief review of the 
MSSM in section 4.2, we present in section 4.3 the allowed region of the SUSY parameter 
space. Section 4.4 shows the various supersymmetric contribution to the b ^ s t^t~ transition, 
while in section 4.5 we study the B decay in some specific SUSY-models, such 

as Minimal supergravity (mSUGRA) model. Minimal flavor violating supersymmetric model 
(MEV) and Extended Minimal flavor violating supersymmetric model (EMEV). As probes of 
new physics effects in 5 —, we propose in section 4.6 to study the ratios Ro{s) 
and R-{s) (introduced and calculated in the SM in Chapter 3) in the generic mSUGRA model. 
Einally, we summarized our analysis in section 4.7. 


4.2 The Minimal Supersymmetric Standard Model 

The minimal supersymmetric standard model (MSSM) is the by far most widely studied poten¬ 
tially realistic SUSY model. It owes its popularity mostly to its simplicity, being essentially a 
straightforward supersymmetrization of the Standard Model (SM), where one introduces only 
those couplings and fields that are necessary for consistency [130-132]. 

The single-particle states of the MSSM fall naturally into irreducible representations of the 
corresponding algebra which are called supermultiplets. Each supermultiplet contains both 
fermion and boson states with the same electric charge, weak isospin, and color degrees of free¬ 
dom, which are commonly known as superpartners of each other. All of the Standard Model 
fermions (the known quarks and leptons) are members of chiral supermultiplets. The names for 
the spin-0 partners of the quarks and leptons are constructed by prepending an “s”, which is short 
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Table 4.1: Chiral supermultiplets in the Minimal Supersymmetric Standard Model. 


for scalar. It seems clear that the Higgs scalar boson must reside in a chiral supermultiplet, since 
it has spin 0. Actually, it turns out that one chiral supermultiplet is not enough. One way to 
see this is to note that if there were only one Higgs chiral supermultiplet, the electroweak gauge 
symmetry would suffer a triangle gauge anomaly, and would be inconsistent as a quantum theory. 

The basic structure of the MSSM is well-known and has been thoroughly discussed in the 
literature [133,134]. We therefore recall just those aspects of the theory which are pertinent 
to 6 —> s transitions. We first display nomenclature conventions for matter superfields 
and their left handed fermion and scalar components in Table 4.1, classified according to their 
transformation properties under the SM gauge group SU (3)c x SU{2)^x1/ (1)y. 

The vector bosons of the SM clearly must reside in gauge supermultiplets. Their fermionic 
superpartners are generically referred to as gauginos. The SU (Sjc color gauge interactions of 
QCD are mediated by the gluon, whose spin-1/2 color-octet supersymmetric partner is the gluino. 
As usual, a tilde is used to denote the supersymmetric partner of a SM state, so the symbols for the 
gluon and gluino are g and g respectively. The electroweak gauge symmetry SU{2)l x U{1)y 
has associated with it spin-1 gauge bosons kk°, W~ and B^, with spin-1/2 superpartners 
iy+, W^, W~ and 5°, called winos and bino. After electroweak symmetry breaking, the kk°, 5° 
gauge eigenstates mix to give mass eigenstates and 7. The corresponding gaugino mixtures 
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Names 

Eermions 

Vectors 

gluino, gluon 

9 

9 

winos, W bosons 

0 

-H 

0 

-H 

bino,B boson 

B^ 

5O 


Table 4.2: Gauge supermultiplets in the Minimal Supersymmetric Standard Model. 


of IU° and are called zino (X°) and photino (7); if supersymmetry were unbroken, they would 
be mass eigenstates with masses mz and 0. Table 4.2 summarizes the gauge supermultiplets of 
a minimal supersymmetric extension of the SM. 

After this brief introduction on the field content of the MSSM, we will review in what fol¬ 
lows, the aspects of the theory relevant to the 6 —> s transitions considering only the case 
of unbroken R-parity. The superpotential which determines the supersymmetry preserving inter¬ 
actions among matter fields is: 


w = psusy HiH2 + YiiQ,um2 + 


(4.1) 


where Q, 77^, L and are the superfields corresponding to the SU (2) doublets and singlets 
for quarks and leptons. Hi and H 2 are the two Higgs superfields, Y^ are the Yukawa matrices and 
/i is the Higgs quadratic coupling. After vector superfield terms are included, the supersymmetric 
Lagrangian schematically appears in component form as [135] 


-^SUSY — asA®- f 


dW 

d^i 


dW 

d^i 


+ 


d^W 

(9$, (9$ 


■iljJCiljj + h.c. 




-V2gG + h.c.] - 


(4.2) 


The index G labels the color, weak isospin and hypercharge factors in the Standard Model gauge 
group, and indices A and B range over the nonabelian subgroups’ adjoint representations. All 
MSSM scalars are assembled into 0, while matter fermions and gauginos are respectively con¬ 
tained within the four-component left handed ifj and A fields. 

Since supersymmetry is manifestly violated in the low energy world, the MSSM Lagrangian 
is supplemented with the soft supersymmetry breaking terms [135] 


C 


soft — 


m-gg^^Cg'^ + + tUsB^CB + h.c. 


— — A^\h\h2 
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- af - e't(Mj)ye5 


+ 


A 1 jqiu'^-h 2 + Af-qi(^hi + + Biisusyhih 2 + h.c. 




(4.3) 


where (mw. and ms), Af {i = 1, 2) and Bfisusy are mass terms for gluino (wino and bino) 
and for the Higgs fields, respeetively. The sealar mass terms M|, M|c, Mj^, M| and M|c are 
in general hermitian 3x3 matriees, while A'^ h 2 , A^ hi and A’^ hi are general 3x3 matriees. 
Allowing all the parameters in (4.3) to be eomplex, we end up with 124 masses, phases and 
mixing angles as free parameters of the model. 

Eleetroweak symmetry breaking induees mixing among MSSM fields. In the matter seetor, 
primed mass eigenstates are related to unprimed gauge eigenstate eounterparts as follows: 


u' = + 

d' = 


U 





d' = 


S^^d \ 
-S^^d-*) 


u' = 
e' = 


z>' = 



(4.4) 


The unitary S and T transformations rotate fermion and sfermion mass matriees into real and 
diagonal forms. The 3x3 quark and lepton mass matriees are simply related to the Yukawa 
eouplings in the superpotential: 


Mu 

Md 

Me 


V2 


V cos (3 






(4.5) 


where vi = v cos (3 and V 2 = v sin (3 are the expeetation values of the two Higgs doublets and 
(/ = Ur,Ul, Dr, Du, Er, Er). The 6x6 squared mass matriees for the squarks (we do 
not present the lepton ones sinee they are not relevant for the fortheoming analysis of SUSY 
eontributions) read: 


M? = T 


u 


+ M2 + ... 

Mu cot (3- 


V2 


fisusyMu cot 13 - 

V2 

+ M2 + ... 




y 


T 


ut 
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Mj = r 

a 


D 


V 


M? + Ml+ ... ^^susyMD tan/3 - 

V2 

tan^( - S‘’’> Af|/5°»' + AfJ + ... 


where M| and M| are the diagonal mass matrices of the up and down squarks and the dots 
stand for terms proportional to m%. Moreover, M| and M| are matrices in the basis in which 
the squark fields undergo the same rotations as the quark ones. This means that we diagonalize 
the matrices Mu and Mo applying rigid rotations to the quark superfields and that there is not 
any flavour change in vertices with both quark and squarks. In the literature this basis is usually 
referred as the SuperCKM (SCKM) one. 

Mixing also takes place in the gaugino and Higgs sectors. The physical Dirac chargino and 
Majorana neutralino eigenstates are respectively linear combinations of left handed Winos, Binos 
and Higgsinos. Thus in the weak eigenstates basis, the chargino and neutalino mass matrices are 
given respectively, by 


and 


M^o 


= Af* 


( 


M^± =W 


rrii 


V 


V2mw sin (3 
V2mw cos [3 -Psusy 


\ 


vt, 


(4.7) 


0 


/ rrig 

0 rriw 

—mz sin 9 cos f3 mz cos 9 cos f3 
\ mz sin 9 sin f3 —mz cos 9 sin f3 


—mz sin 9 cos f3 
mz cos 9 cos f3 
0 

fJ'susy 


mz sin 9 sin f3 \ 
-mz cos 6* sin f3 


fJ'susy 

0 


J 


(4.8) 


where, the unitary transformations U, V and M are unitary matrices, which diagonalize these 
fields mass matrices. 

After the gauge eigenstate fields in the supersymmetric Lagrangian 4.2 are rewritten in terms 
of their mass eigenstate counterparts, * it is straightforward to work out the interactions of 
gluinos, charginos and neutralinos with quarks and squarks. We list below the resulting terms 
which participate at one-loop order in di dji+i- decay: 


8 

(4.9) 

a=l 

2 4 

+ + Zl^^R)d + h.c., (4.10) 

i=i i=i 


We suppress primes on mass eigenstate fields from here on. 



79 


where 




2;Dl 


= 92 




Mr, 


72 - 


= 92 


Unr^^V, 


\/2mw sin /9 
Mo 


Vc 


CKM, 


CKM 


\/2mw cos (3' ’ 


92 

V2 


(-W 2 + 5tariW;,)T“‘ +W3r°« „ 

3 TTZiyCOSp. 




Mr 


, I ^/- (4.11) 

-\/2 Vo mw cos py 

We see clearly that the flavor mixing enters into these interactions through the Kobayashi- 
Maskawa matrix Vckm = and the 6x3 block components of and T°: 


u{D) _ f^UUDo) rj-UR{Du) 
-*6x6 ~ I -*6x3 ! -*6x3 


(4.12) 


The Feynman rules for all these interactions may be found in the literature [133,136]. Having set 
up the basic MSSM framework, we are now ready to explore its large parameter space. We take 
up this topic in the following section. 


4.3 MSSM parameter space 

Before predictions can be derived from the Minimal Supersymmetric Standard Model, explicit 
values for the parameters in the superpotential (4.1) and soft supersymmetry breaking Lagrangian 
(4.3) must be specified. 

In order to determine the allowed region of the SUSY parameter space, we require the fol¬ 
lowing phenomenological constraints [36] 

(1) 6 ^ S 7 constraint from BABAR [11], CLEO [12], ALEPH [13] and BEEEE [14] collab¬ 
orations, i.e., 3.06 X 10-^ < B{B X^p) < 3.85 x 10“^ 

(2) Erom the recent experiment at EEP 2 [36], we impose that all the charged SUSY particles 
are heavier than 70 GeV 

(3) The neutarlino x'j’ mass is larger than 32 GeV. 

(4) All sneutrino masses are larger than 43 GeV [36]. 

(5) The gluino and squark mass bounds from Eermilab Tevatron experiments [137]. The pre¬ 
cise bounds on the gluino mass and the averaged squark mass except for the top squark is 
restricted to be larger than about 180 GeV. 

(6) The stop ii mass is larger than 87 GeV. 

Having these phenomenological constraints at hand, let’s explore the SUSY contributions to 
the b —s- transition. 
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4.4 Supersymmetric contribution to s £+£- 

The Major theoretical breakthrough in the analysis of the ECNC tests in SUSY models came in 
1983 when Duncan [139] and, independentely, Donoghue, Nilles and Wyler [140] noticed that 
EC transitions could occur at the g — q — q vertices. The replacement of the weak coupling of 
the W with the strong coupling of the gluino, and the presence of approximately the same CKM 
mixing angles, raised the hope of possible SUSY enhancements of ECNC processes through one- 
loop gluino exchanges. Indeed, in the general context of MSSMs, it was found that these gluino 
mediated ECNC contributions could play a relevant role in CP violation of the K-system [141] 
and, more recently, in the radiative b ^ sy [142] and b ^ sg [143] decays. 

Within the MSSM, the b ^ s transition is governed by five possible classes of contri¬ 
butions. They correspond to four classes of one loop SUSY diagrams that produce b si+i- 
transition, in additions to W-exchange. We will classify them according to the particles running 
in the loop: 

(1) W~ and up quarks (SM), 

(2) H~ and up quarks, 

(3) and up squarks, 

(4) g and down squarks, 

(5) and down squarks. 


We list below the VE-scale matching contributions to the Wilson coefficients Cj, Cg and 
Cio which arise from one-loop MSSM diagrams, presented in Pig. (4.1) and (4.2). The total 
contribution to Ci {i = 7, 9,10) and the normalized ratio can thus be written respectively as 

Cf\mw) = C^^{mw) + Cr{mw), (4.13) 


and 


Cf^imwY 


(4.14) 


We start with the W -scale matching contributions to the coefficient C 7 of the magnetic mo¬ 
ment operator in the AB = 1 effective Hamiltonian which arise from one-loop MSSM diagrams, 
presented in Pig. (4.1) (only the 7 -penguin diagrams). Thus its total contribution reads as 


C\°^{mw) = C^^irnw)+ CY^{rnw), 


(4.15) 


where CY^{mw) represent the new physics contribution at the scale mw to the U 7 , defined as: 

CY^{mw) = + bcf {niw) + 5Cf {niw) + 5Cl{mw)- (4.16) 
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The explicit expressions for the various terms are [135]: 
• Standard Model graphs: 


c-f" = ^Mx.) 


Graphs with charged Higgs loops: 


sc!;* =iii^coei3M^ 


6 1 2 m 


/!=•= 


/!=•= ’ 


+ / 2 (: 


\± ■ 


(4.17) 


(4.18) 


Graphs with chargino loops: 
1 


6C^ = 


6 2 2 


^glv^iVtb “ “ 


2 


(4.19) 




m. 


UA 


2 

Xi 




rrife 


2 

X/ 


Graphs with neutralino loops: 

^ 6 4 


= 


EE 


m 


w 


^92Vt*sVtb ^ ^ rri-o 


(4.20) 




m 




m-0 

X? 


0 4. n 

+ (2F*)L(2F")A.^f4 


mb 


2 

TTlj 

m?o 

X? 


Graphs with gluino loops: 


SCI = 




6 2 
^9 


^9iV,lVtb ^ mj 

yC(T*'*)UT*'*: 


(4.21) 


A3 


/s 


m 


dA 


m~ 


+ 2{T*’4)l{T°‘) 


™3 f 

A3 J4 

mb 


m 


dA 


ml 


The one-loop integral functions which enter into these matching conditions are given in Ap¬ 
pendix C. Concerning the SUSY contributions to the chromo-magnetic coefficient Cg, it has the 
same structure as the Cj ones, with different colour factors and loop-functions. The lU-scale 
matching contributions to the semileptonic coefficients , namely and Gio, which arise from 
one-loop MSSM diagrams, presented in Figs. (4.1) and (4.2), and can be written respectively as 


^ 9^0 i^w ) — i^w ) + CgJIf {mw ), 

where the SUSY contributuions to the semileptonic coefficients, are defined as: 


(4.22) 


iZ,'y)-peng,(H,x,X°,9) , xn'Cn)-peng,{H,x,C ,9) , xm(x,X°)-box xm'C^O-box 


C^^^imw) = +5Cg 


— box 


(4.23) 
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Eigure 4.1: SUSY-Penguin diagrams relevant to the semileptonic b ^ s transition in the 
MSSM. The cross denotes a possibility to attach the photon or the Z-boson. 


The various terms are given as follow [135]: 

• Z-penguin graphs with eharged Higgs loops: 


^ ^ 3 ,oei3x,M^) 

8 sm Ow 

^C'Z-peng,H ^ 1 ^^^2 


9 n , 

cot^ f3xt hi— 
2 m 


'h± 


7 -penguin graphs with eharged Higgs loops: 

j:^''y-peng,H _ 1 

“ 18 


'h± 


m 


Z-penguin graphs with ehargino loops: 

c^Z-peng,x _i 

OWn in — =■= 


X 


2t/| 




A,B=1 I,J=1 


(4.24) 


(4.25) 
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Figure 4.2: SUSY-Box diagrams relevant to the semileptonic b ^ s transition in the MSSM. 
The bubble in figure {d) reflect the Majorana nature of the neutralinos. 


+ 2 ^xf ^Xj ^0 Ida U 


Xi’ Xj' 


92^tsytb ^_ b = 17 , j =1 


x|c2(m2±,m?^,m?^) 5/j - C2{ml^,m^^±,wf±) 5as Vn 


+ \rn^± m-± Co(m?^,m^±,m^±) SabUjiU 


Xi ' Xj‘ 


(4.26) 


7 -penguin graphs with ehargino loops: 


= 0 


^QX-peng,x 


6 2 n 

mr, 






(4.27) 


Z-penguin graphs with neutralino loops: 


6 4 


SCZ-peng,x° = - V V (2f^)L(^ 


J )Bi 


X |c2(m|o, )as ()/j - C2(m|^, r«|o, m|o) ()as (A/'/sA/'js - 

-i m^o m^o co(m|^,m|o,m|o) ()as (MsA/j'g - A/)4A/';4)| 
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^(j'Z-peng,x° _ 




2A\^J JBi 


^2 ts a,b=1i,j=1 

X {c2(m|o, Sjj - C 2 (m|^, m|o, m|o) 5ab (Ks-^js - K4-^j4) 

-i m^o m^o Co (m|^, m|o, m|o) S^b Wiz-^jz “ | (4.28) 

7 -penguin graphs with neutralino loops: 
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(4.29) 


Z-penguin graphs with gluino loops: 
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7 -penguin graphs with gluino loops: 
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Chargino box graph: 
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Neutralino box graphs: 
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The one-loop integral functions which appear within these MSSM matching conditions are given 
in Appendix C. The total Wilson coefficients given above, namely and evaluated 

perturbatively at the W scale and then evaluated down to the renormalization scale n h by 
the renormalization group equation (RGE). The details of this strong interaction running are quit 
cumbersome, and we will not present them here. However, the details of this RGE can be found 
in [20], 


4.5 Analysis in supersymmetry 

Having now at hand the general MSSM contribution to the b ^ s transition, we will turn 
now to more restricted framework of the MSSM. We employ the following models to study the 
rare B decays: 

(1) Minimal supergravity (mSUGRA) [144-146], 

(2) Minimal flavor violating supersymmetric model (MEV) [147], 

(3) Extended Minimal flavor violating supersymmetric model (EMEV) [148]. 

The last of these models serves as a generic supersymmetric extension of the SM having non- 
CKM flavor violations. 


4.5.1 SUGRA model 

In supergravity (SUGRA) model [149,150], which is a result of the unification of the super- 
symmetry transformations with the space-time symmetries of general relativity, the soft SUSY 
breaking terms are supposed to arise from a hidden sector of the theory which can only commu¬ 
nicate with the ordinary matter fields through gravitational interactions. 

Since gravity is flavour-blind, the breaking terms can be realized in a minimal version [151], 
introducing a common scalar mass parameter mg and trilinear coupling Ax, a universal gaugino 
mass parameter M 1/2 and the bilinear Higgs parameter b. This reduced set of breaking parameters 
is called minimal supergravity (mSUGRA). 

In the minimal SUGRA model, the soft SUSY breaking terms are assumed to take the fol¬ 
lowing universal structures at the Planck or GUT scale: 


(M|)l 

A? 


= (MDl = {Mf)l = {MDI = , 

A^ = 

flAxmo , = fiAxmo , = fjAxmo , 

niw = mg = Mgx . 


(4.34) 
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In the minimal supergravity model the soft breaking parameters mo and Aq are assumed to be 
equal^ With the above initial conditions we can solve the one-loop RGBs for the SUSY breaking 
parameters and determine these parameters at the electroweak scale [152]. We also require that 
the electroweak symmetry breaking occurs properly to give the correct boson mass. 

Scanning over the soft SUSY breaking parameter space in the range mo < 600 GeU, Aq < 
600 GeU, Mgx < 600 GeU and \Ax\ < 5 for each fixed value of tan{3, it turns out that the 
parameter space of this model may be decomposed into two qualitatively different regions: [31, 
145] 


• Eor small tan (3, say tan /9 ~ 3, the sign of is the same as in the SM. Hence, no spec¬ 
tacular deviations from the SM can be expected in the B K* f+r decay mode. Given 
the theoretical uncertainties in the SM estimates, it would be very difficult to disentangle 
any SUSY effects for this scenario in this decay [31]. 

• Eor large tan f3, the situation is more interesting due to correlations involving the branching 

ratio for B the mass of the lightest CP-even Higgs boson, ruh, and sign(/is^, 5 y). 

In this case, there are two branches for the solutions for rrih and B{B Xs'j). The 
interesting scenario for SUSY searches in B ^ is the one in which sign{jj,susy) 

and TUh admit to be positive^. 

However, a large^ tan (3 values can get a considerable impact on the 6 — s branching ratio: 

• At the low-s region the branching ratio B —> fj,~ could be enhanced by about 30% 

compared to the SM one, nevertheless this enhancement is difficult to disentangle from the 
non-perturbative uncertainties attendant with the SM-distributions [31]. 

• The B dilepton mass distribution could be enhanced by about 100% and 

this is distinguishable from the SM-related theoretical uncertainties [31]. A very similar 
supersymmetric effects have been worked out for the inclusive decays B [145], 

where enhancements of (50-100)% were predicted in the low-s branching ratios. 

Summarizing for the SUGRA theories, small tan [3 implies the sign of G 7 ®® being the same as 
in the SM. Hence, no spectacular deviations from the SM can be expected in the B K* f+r 
decay mode. However, large tan 13 solutions lead to G 7 ®® being positive, and one expects an 
enhancement up to a factor two in the dilepton mass distribution in B ^ This would 

be a drastic deviation from the SM, which cannot be fudged away due to non-perturbative effects. 

Whereas in the nonminimal case we treat the two as independent parameters 

Wor example, this happens for tan (3 > 10, in which case nih = (115-125) GeV and Gy®® is positive and obeys 
the B Xgj bounds [146]. Following the generic case shown earlier, one expects a constructive interference of 

the terms depending on Gy®® and Gg in the dilepton invariant mass spectra. 

Han/3 = 30 corresponds to [146]: ry =—1.2, rg = 1.03, rig = 1.0 . 
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Figure 4.3: 90% C.L. bounds in the [R^{|Ji), Rs^n)] plane following from the world average 
B —> Xs 7 branching ratio for p = 'rny/ (left-hand plot) and p = 2.5 GeV (right-hand plot). 
Theoretical uncertainties are taken into account. The solid and dashed lines correspond to the 
rUc = me,pole and me = m^^(ph) cases respectively. The scatter points correspond to the 
expectation in MFV models (the ranges of the SUSY parameters are specified in the text) [20]. 


4.5.2 The MFV model 

The minimal flavour violating (MFV) SUSY model [147] is based on the assumption of minimal 
flavor violation, whieh means that all the genuine new sourees of flavour ehanging transitions 
other than the CKM matrix are switehed off. Here, quarks and squarks are aligned so there is 
no flavor-ehanging q — q' — (Z, 7, g) vertex and the eharged one, d — u — is governed by 
the CKM matrix. As a eonsequenee, in this model neutralino-down-squark and gluino-down- 
squark graphs do not eontribute to either 6 —> 57 or 6 —^ transitions. In addition to the 

charged Higgs-top graphs, chargino-up type squarks loops with a light stop U, and the lU^-top 
quark loops, present in the SM, give the dominant contribution. While not holding generally, the 
assumptions in the MFV-SUSY model are valid over an important part of the MSSM parameter 
space. 

In order to derive the new physics range contributions to Cg and Cig, we scan over the pa¬ 
rameter space in the range 78.6 GeU < mH± < 1 TeU, 0 < \psusy\ < 1 TeU and 
THq = 771^2 ^ 90 GeU, where denotes the (degenerate) masses of other than top squarks, and 
mf, > 50 GeU. We reject too light charginos, demanding m ± > 90 GeU and also taking the 
B —> Ys 7 experimental constraint. We have chosen a stop mixing angle < 7r/2, i.e. the light 
stop [2 = —sinOdL + cosOdji is almost left handed. In order to produce bounds that can be 
compared with the model independent allowed regions plotted in Fig (4.4), it turns out that the 
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surviving SUSY points can be divided to two sets [20]: 

• Eor small tan/?, for which we again take tan/? = 2.3, we find that the ratio i ?7 remains 
positive, i.e. < 0, and lies within the experimentally allowed bounds from B Xs'y, 
and the corresponding bounds on the semileptonic coefficients are in the range — 0.2 < 

< 0.4 and -1.0 < < 0.7. 

• Eor large tan [3, taken to be 50, C 7 ®® changes sign (R^ < 0). The corresponding bounds 

on and tend to be in the range —0.2 < {^w) < 0.3 and —0.8 < 

< 0.5. 


The above discussion applies to any supersymmetric model with flavour universal soft-breaking 
terms, such as mSUGRA MSSM and gauge-mediated supersymmetry breaking models. Beyond- 
the-SM flavour violations in such models are induced only via renormalization group running, 
and are tiny. Before ending this subsection, let us discuss the impact of 6 —57 on MEV models 
in varying the MEV SUSY parameters: 

• The strong correlation between the values of the Wilson coefficients C 7 and Cg, shown in 
Eig. 4.3. In fact, the SUSY contributions to the magnetic and chromo-magnetic coeffi¬ 
cients differ only because of colour factors and loop-functions. 

• The dependence of the charged Higgs contribution to magnetic coefficient at the scale Hb, 
can be seen from Eig. 12 in ref [20]. It turns out that with a specific scenario, it is possible 
to obtain lower bounds on some SUSY particles. 

• The chargino contributions on Ci{iJib), show very strong consequences^ on the b ^ S'j 

transition as it has been shown in [20]. Assuming for instance rrif^ = = 500 GeV 

the chargino contributions (normalized to (sin 6 *^ tan (3)) is of order 0.2. If one then allows 
for larger values of the stop mixing angle and of tan (3, the contribution can easily violate 
the b ^ sj constraint by more than one order of magnitude (e.g. for sin 6 *f = 0.5 and 
tan (3 = 50 one obtains something of order 6 that is orders of magnitude above the current 
limit) [ 20 ]. 

Up to this point we have been looking at the minimal case where all genuine sources of flavour 
changing transitions in the MSSM are attributed to the CKM matrix elements. However, an 
interesting question to ask is whether there might be other, non-universal scenarios, where new 
flavour changing transitions (other than the CKM mixing elements) could occur. A model which 
incorporates these features is the so-called Extended-MEV (EMEV) model [148], discussed in 
the following subsection. 

^As it has been argued in [20], one can exploit the 6*7 and tan (3 dependence since (for non negligible values of 
the stop mixing angle) the chargino contribution is essentially proportional to sin 0^ tan (3. 
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Figure 4.4: NNLO Case. Superposition of all the constraints. The plots correspond to the 
Af^{2.5 GeV) < 0 and A^f^(2.5 GeV) > 0 case, respectively. The points are obtained by 
means of a scanning over the EMFV parameter space and requiring the experimental bound 
from B —> Xgy to be satisfied [20]. 


4.5.3 The Extended-MFV model 

The Extended-MFV (EMFV) model whieh is a generalization of the MFV-model using the lan¬ 
guage of minimal insertion approximation (MIA) [153] in a supersymmetrie context. Its main 
assumption lies in the fact that new sources of flavour changing transitions other than the CKM 
mixing elements could occur. 

EMFV-models are based on the heavy squarks and gluino assumption. Moreover, the charged 
Higgs and the lightest chargino and stop masses are required to be heavier than 100 GeV in order 
to satisfy the lower bounds from direct searches. The rest of the SUSY spectrum is assumed to be 
almost degenerate and heavier than 1 TeV. The lightest stop is almost right-handed and the stop 
mixing angle (which parameterizes the amount of the left-handed stop ii present in the lighter 
mass eigenstate) turns out to be of order ~ 10%; for definiteness we will take \9i\ < vr/lO. 

The assumption of a heavy (> 1 TeV) gluino totally suppresses any possible gluino-mediated 
SUSY contribution to low energy observables. Note that even in the presence of a light gluino 
(i.e. rrig ~ 0(300 GeV)) these penguin diagrams remain suppressed due to the heavy down 
squarks present in the loop. On the other hand, the presence of only a single light squark 
mass eigenstate (out of twelve) has strong consequences due to the rich flavour structure which 
emerges from the squark mass matrices. Adopting the MIA-framework [153], all FC transitions 
which are not generated by the CKM mixing matrix are proportional to the properly normalized 
off-diagonal elements of the squark mass matrices: 




(4.35) 
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where j = 1,2,3 and A,B = L, R. In this approach, some remarks are in order: 


• The only sizable contributions arise from the inserted mass insertions involving the light 
stop. 

• All the other diagrams require necessarily a loop with at least two heavy (> 1 TeV) squarks 
and are therefore automatically suppressed. 


This leaves us with only two unsuppressed flavour changing sources other than the CKM matrix, 
namely the mixings ul — h (denoted by and cl — h (denoted by We note that 
and are mass insertions extracted from the up-squarks mass matrix after the diagonalization 
of the stop system and are therefore linear combinations of (5i3)i/j, ('^i3)ll i^23)LR^ 

{^ 23 )ll^ respectively. The insertions relevant to our discussion are normalized as follows: 


5 , 


u{c)i,i2 


^ ^u{c)Li2 \ Vtdis)\ 

K*d{s) 


(4.36) 


The phenomenological impact of has been studied in [148] and its impact on the b ^ s'y 
and b si^i~ transitions is indeed negligible. Therefore, we are left with the MIA parameter 
Thus, the SUSY parameter space that we have to deal with are the same as the MEV- 
model with MIA parameter In order to explore the region in the [Cg^, plane (where 
are the sum of MEV and MI contributions) that is accessible to these models, one has to 
apply the b ^ s'j constraint on the EMEV parameter space, which are the same as the MEV-ones 
apart from \6i\ < tt/IO and The surviving points are shown in Eig. 4.4 together with 

the model independent constraints. To get them, one has to use the integrated branching ratios to 
put constraints on the effective coefficients. This procedure allows multiple solutions, which can 
be disentangled from each other only with the help of both the dilepton mass spectrum and the 
forward-backward asymmetry. Only such measurements would allow us to determine the exact 
values and signs of the Wilson coefficients Cj, Cg and Uiq. 


4.6 The ratios R-{s) and i?o(s) 3S probes of New Physics in 

B K *t+t 

In order to look for new physics m B ^ we propose to study the ratios Ro{s) and 

R-{s), introduced in the previous section. As well known, new physics can distort the dilepton 
invariant mass spectrum and the forward-backward asymmetry in a non-trivial way. 

To illustrate generic SUSY effects in i? —we note that the Wilson coefficients 
U|^, Cg and Uio receive additional contributions from the supersymmetric particles. We 
incorporate these effects by assuming that the ratios of the Wilson coefficients in these theories 
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Figure4.5: The Ratios R -{ s ) (left-hand plot) and Ro { s ) (right-hand plot) in the Standard Model 
with Rb = 0.094, Csc/(3) = 1-3, ^1^‘Ho) = 0.28 and in SUGRA, with ( ry , rg) = (1.1, 1.4). 
The SM and the SUGRA contriburions are represented respectively by the shaded area and the 
solid curve. The shaded area depicts the theoretical uncertainty on Csu{ 3 ) = 1-3 ± 0.06 and on 
= 0.28 ± 0.04 [29]. 

and the SM deviate from 1. These ratios for k = 7,8, 9,10 are defined as follows^: 

nSUSY 

= -^mT • ( 4 - 37 ) 

They depend on the renormalization seale (exeept for Ciq), for whieh we take p = nib^poie- For 
the sake of illustration, we use representative values for the large-tan /3 SUGRA model, in whieh 
the ratios rj and rg aetually keep or ehange their signs. The supersymmetrie effeets on the other 
two Wilson eoeffieients Cg and Gio are generally small in the SUGRA models, leaving rg and rio 
praetieally unehanged from their SM value. To be speeifie, we take the two allowed seenarios* ** 

ry = -1.2, rg = -1, rg = 1.03, rig = 1.0, (4.38) 

and 

ry = 1.1, rg = 1.4, rg = 1.03, ryg = 1.0. (4.39) 

In Figs. (3.27), (3.28), (4.5) and (4.6), we present a eomparative study of the SM and SUGRA 
partial distribution for iT_ and Hg, respeetively. In doing this, we also show the attendant theo- 
retieal uneertainties for the SM, worked out in the LEFT approaeh [28]. For these distributions, 

II These ratios have been introduced in Eq. (4.14), and we change slightly their notations by instead of Ri and 
keep this later notation for our helicity ratios. 

** We thank Enrico Lunghi for providing us with these numbers. 
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Eigure 4.6: The Ratios (s) (left-hand plot) and Rq{s) (right-hand plot) in the Standard Model 
with Rb = 0.094, Csc/(3) = 1-3, ^1^‘^(0) = 0.28 and in SUGRA, with {ry, rg) = (-1.2, -1). 
The SM and the SUGRA contributions are represented respectively by the shaded area and the 
solid curve. The shaded area depicts the theoretical uncertainty on Csu{3) = 1-3 ± 0.06 and on 
= 0.28 ± 0.04 [28]. 


we have used the form faetors from [31] with the SU(3)-symmetry breaking parameter taken in 
the range Csu{3) = 1-3 ± 0.06. 

Erom Eig. (4.5) , where >0, (k = 7, 8, 9,10), it is diffieult to work out a signal of new 
physies from the SM pieture. There is no surprise to be expeeted, due to the faet that in these 
seenario the corresponding ratio is approximatively one, which makes the SUGRA picture 
closer to the SM one. However, Eig. (4.6) with (rj, rg) < 0 illustrate clearly that despite non- 
perturbative uncertainties, it is possible, in principle, in the low s region to distinguish between 
the SM and a SUGRA-type models, provided the ratios differ sufficiently from 1. 


4.7 Summary and Outlook 

In this chapter we have presented a phenomenological profile of the semileptonic rare B —> 
decay in the context of supersymmetric theories. Considering the straightforward su- 
persymmetrization version of the SM, the so-called MSSM, we have reviewed in details its con¬ 
tribution to the b process. We illustrate the constraints on the Wilson coefficients C 7 , 

Cg, Cq and Cio appearing in the effective Hamiltonian formalism, that the current data on rare B 
decays implies in the context of minimal flavour violating model and in more general scenarios 
admitting additional flavour changing mechanisms. Einally, incorporating these supersymmetric 
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effects on the corresponding Wilson coefficients, we have shown their phenomenological impact 
on the ratios Roi-s) and i?_ (s) as probes of new physics. 

Our studies can be reported as follows [28] : 

• We have shown within the MSSM framework the complete SUSY contributions to the 
b si+i- decay. Beyond the lU-exchange in the SM, four other classes contribute to 
this process, namely the Higgs-exchange, Chargino-exchange, gluino-exchange and the 
Neutralino-ones. 

• Using the current data constraints on rare B decays, we have discussed their phenomeno¬ 
logical impacts on the Wilson coefficients U7, Ug, Ug and Uiq. It turns out that present 
experimental measurements leave considerable room for beyond-the-SM contributions , 
especially, to the allowed region in the [Cg, Uio] plane. However, in the SUGRA models 
they are practically unchanged from their SM value. 

• In order to look for new physics in B ^ K* f+r ', we propose to study the ratios Ro{s) 
and R_{s). For the sake of illustration, we use representative values for the large-tan/? 
SUGRA model with two scenarios, namely (r7, rg) < 0 and {rj, rg) > 0 . We have 
noticed that despite non-perturbative uncertainties, it is possible when (r7, rg) < 0 , in 
principle, in the low s region to distinguish between the SM and a SUGRA-type models, 
provided the ratios differ sufficiently from 1. 




Chapter 5 


Summary & Future 


While waiting for the completion of the second-generation experiments at hadron colliders, 
BTeV (Fermilab) and LHCb (CERN), 5-physics is among the most active and promising fields 
in recent particle physics. Its importance lies in the deeper understanding of the Standard Model 
and particle physics in general. The futur investigation of b decays (such as rare 5-decays, two 
leptonic 5-decays, ■ ■ ■) at the 5-factories, namely BABAR and BELLE, and at hadron colliders 
will probe the flavour sector of the SM with unprecedented precision, and may be reveal new 
physics effects. 

Rare 5 decays involving flavour-changing-neutral-current (ECNC) transitions, such as 6 —> 
sy and b , have received a lot of theoretical interest [154]. Especially, after the first 

measurements of the radiative decay 5 ^ Xsy were reported by the CLEO collaboration [155] 
in 1995 and recently the first measurements of the semileptonic rare 5-decays reported in the 
inclusive (exclusive) 5 — Xs{K, by the BELLE collaboration [14] (BABAR collabo¬ 

ration [18]). With increased statistical power of experiments at the 5-factories in the next several 
years, the decays discussed above and related rare 5 decays will be measured very precisely. We 
summarize the projections for improvement in the experimental contributions to the precision 
of CKM matrix elements Vuh, Vch and Vts [156] in Table 5.1 while in Table 5.2 [156] the decay 
reach of the rare 5-decays is given. 

Within this thesis, we have reported an 0(Q;5)-improved analysis of the various helicity com¬ 
ponents in the decays 5 —and 5 — piu^, carried out in the context of the Large- 
Energy-Effective-Theory. Our studies can be summarized as follows [28,29] : 

• The underlying symmetries in the large energy limit lead to an enormous simplification 
as they reduce the number of independent form factors in these decays. The LEET- 
symmetries are broken by QCD corrections, and we have calculated the helicity com¬ 
ponents implementing the 0{as) corrections. 

• The results presented here make use of the form factors calculated in the QCD sum rule 
approach. The LEET form factor ' (0) is constrained by current data on 5 ^ K*y. 
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Experimental 

Measurement 

a{%) 

2001 

stat/sys 

a(%) 

2006 

stat/sys 

a(%) 

2011 

stat/sys 

Vub 

BiB ^ p iue) 

4.3/8 

8.6/2.4 

1.4/2.4 


B{B —u jvf) 

3.4/16 

4.0/2.4 

2.8/2.4 

Vcb 

B{B D iue) 

3.1/4 

0.4/2 

0.10/1 


B{B c iue) 

2.5/2 

0.3/1 

0.07/0.5 

Vts 

Arris 





Table 5.1: Projections for improvement in the experimental contributions to the precision ofCKM 
matrix elements Vuh, Vcb and Vts [156], 

• As the theoretical analysis is restricted to the lower part of the dilepton invariant mass 

region in B ^ typically s < 8 GeV^, errors in this form factor are not expected 

to severely limit theoretical precision. This implies that distributions involving the H_ (s) 
helicity component can be calculated reliably. Precise measurements of the two LEET 
form factors (s) and (s) in the decays B —> can be used to largely reduce the 
residual model dependence. 

• With the assumed form factors, we have worked out a number of single and double (Dalitz) 
distributions in B ^ piui, which need to be confronted with data. 

• An analysis of the decays B is also carried out in the so-called transversity 

basis. We have compared the EEET-based amplitudes in this basis with the data currently 
available on B ^ K* J/ip and find that the short-distance based transversity 
amplitudes are very similar to their long-distance counterparts. 

• We also show the 0{as) effects on the forward-backward asymmetry, confirming essen¬ 
tially the earlier work of Beneke, Feldmann and Seidel [25]. 

• Combining the analysis of the decay modes B and B —> piu^, we show that the 

ratios of differential decay rates involving definite helicity states, R-{s) and Ro{s), can be 
used for testing the SM precisely. We work out the dependence of these ratios on the CKM 
matrix elements |Eu6|/|V[s|. We have also analyzed possible effects on these ratios from 
New Physics contributions, exemplified by representative values for the effective Wilson 
coefficients in SUGRA models. 

The main thrust of this work lies, however, on showing that the currently prevailing theoret¬ 
ical uncertainties on the SM distributions in B ^ can be largely reduced by using the 

EEET approach and data on B ^ K*'^ and B decays. 
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CHAPTER 5. SUMMARY & FUTURE 




Hadron Collider Experiments 

e'''e 5 Factories 

Decay Mode 

Branching 

CDF 

BTeV 

ATLAS 

BaBar 

Super- 


Fractions 

DO 

LHCb 

CMS 

Belle 

BaBar 



(2fb-^) 

(loX) 

(1 Year) 

(0.5 ab“^) 

(10 ab“^) 

5^X,7 

(3.3 ±0.3) 




IIK 

220K 


X 

o 

1 




1.7K 

34K 






(B Tagged) 

(B Tagged) 

5 ^ K*y 

5 X 10“^^ 

170 

25K 


6K 

120K 

5 ^ p{u)y 

2 X 10“^^ 




300 

6K 

5 ^ Xsp^p~ 

(6.0 ±1.5) X 10“^^ 


3.6K 


300 

6K 

5 ^ X,e+e- 





350 

7K 

5 — > K*p^p~ 

(2± 1) X 10“® 

60-150 

2.2K/4.5K 

665/4.2K 

120 

2.4K 

5 ^ K*e+e- 





150 

3K 


Table 5.2: Decay reach ofB experiments for rare decays [156]. 

Finally, we remark that the current experimental limits onB ^ (X^, K, decays [ 157 - 

160 ] are already probing the SM-sensitivity. With the integrated luminosities over the next couple 
of years at the 5-factories, the helicity analysis in 5 —piug and B decays presented 

here can be carried out experimentally. This work will help the search for flavour changing neu¬ 
tral current B —> and in particular, will contribute to precise determinations of the 

FEET form factors, and the CKM matrix elements |Uub|/|Vts| using as well the B —> piu^ decay 
in forthcoming 5-facilities. 
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Appendix A 


Generalities 


A.l Input Parameters 


M\y 

80.4 GeV 

fs 

200 ± 20 MeV 

rhtimt) 

167 ± 5 GeV 

fKW\ 

225 ± 30 MeV 

^&,po/e(2 GcV) 

4.6 ±0.1 GeV 

fK*,±i^GcY) 

185 ± 10 MeV 

nic 

1.4 ± 0.2 GeV 

fp (1 GeV) 

198 ± 7 MeV 

^em 

1/137 

^B,+ 

(3 ± 1) GeV 

Tb 

1.65 ps 


0.2 ±0.1 

KV,i\ 

0.041 ±0.003 


0.05 ±0.1 


0.094 ±0.014 

±■’(0) 

0.28 ±0.04 

A (^/=5) 

^'■QCD 

220 ± 40 MeV 

±’(0) 

0.22 ±0.04 

{«;■)? 

0.3 ±0.2 (GeV)-i 

{s-'>[r’ 

3.48 


Table A.l: Input parameters and their uncertainties used in the calculations of the 
decay rates for B and B —> piui in the LEET approach. 








99 


A.2 QCD 


The QCD Lagrangian reads in covariant gauge [161] (A": gluon field) 

^^QCD = ^ — ^q)(l — + ^fix + ghosts , 

q=u,d,s,c,b,t 

where the gauge fixing term and the one for ghosts c“, c“ are given as 

The chromomagnetic field strength tensor and covariant derivative are written as 

G% = d,Al-d,Al + gryAlAl, 

= d,-tgT^Al, 

Df = 5^^d,-gr^^Al, 

where are the structure constants of SU (3), defined by 

j'^a _ j^j^abxrjnx 

We have the identities 

[D„D,] = -igrci,, 

D;^Df-D:<’Df = -gr'^Gl,. 

Often the abbreviation is used 


(AT) 


(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 


G,, = . (A. 10) 

T“, a = 1,..., 8 are the generators of QCD. They are related to the Gell-Mann (3x3) matrices 
A“ through T“ = The T“ obey the following relations (f, j, /c = 1, 2, 3) 


SviTa) = 0, (A.ll) 

Sp{TJ',) = 5abl2, (A.12) 

Wu = + , (A,13) 

TSTt, = SitG, (A.14) 

with the invariant Cf in an arbitrary SU{Nc) given as 


The coefficients of the QCD beta function (see Eq. (2.29)) are written as: 

/5o - -^-, 

34Ar2 _ lOAT^AT^ - eCirA^/ 


(A. 16) 
(A. 17) 


Here, N^. denotes the number of colours (A), = 3 for QCD) and Nf denotes the number of active 
flavours {Nf = 5 for the effective Hamiltonian theory relevant for b decays). 
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APPENDIX A. GENERALITIES 


A.3 Feynman Rules 

The CO variant derivative consistent with our definition of the operator basis and the corresponding 
Wilson coefficients given in section 2.2 is [53] 

D^ = d^ + igT^^Al + ieQA^ , (A. 18) 

where A°^^ A — fi denote the polarization four-vectors of the gluon, photon respectively. Note 
that the sign convention of the strong coupling here is opposite to the usual one appearing in 
QCD text books [161,162] given in Eq. (A.5), but can be made consistent with the substitution 
g — —g. The Eeynman rules consistent with eq. (A. 18) are given here with boson propagators 
in Eeynman gauge. In a general gauge with gauge parameter ^ they are written as: 

.g^.u + (^ - (k'^ + ie) , ^ 

with .^ = 1,0 corresponding to Eeynman, Landau gauge, respectively. The Eeynman rules are : 


q 


j k 


^ -S. 




7, .9 

vwwwwv 

a, ji 


h, V 




//• I 7, .9 

— ^ - 

:i k 




complemented by the rules : 

• evaluate fermion lines against the momentum flow 

• add a (— 1 ) for a closed fermion loop and perform the trace over the string of 7 matrices 

The rule for an O 7 operator insertion is, using = iq^ for an out going photon and further 
e ■ g = 0 for a real photon, and = d^A’' — d'' A^^, 

aF = 4] = . (A.20) 

The Eierz transformation in d = 4 dimensions is defined as: 


{(lll^^Lq2){q3l|^Lq4 = +{qi7f,Lq4{q^jf,Lq2) , 

{q^^^L{R)q2){q3^^R{L)q4 = {-2){q-^R{L)q4{q^L{R)q2) . 


(A.21) 

(A.22) 
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A.4 Utilities 

A variety of tools for 1-loop calculations is collected in the appendix of ref. [162]. 
Distributions: 

(5(a;) = — [ dqe^'^^dq , 

27r Jr 

—i f 

9{x) = lim — / dq - dq , 

2n Jr q — ie 


Geometrical series : 


dx 


1 ± a; 


{x) = S{x) 


lTx + x‘^ + y^(=F 


X] 


n=3 


Special Functions useful for loops 
Poly-logarithms : 


Spence function : 


oo u 


Lin{z) = kl < 1, 


k=l 


dt 


Li 2 {z) = -/ ^ln(l-t). 
Jo t 


Sp{z) = Li 2 {z) = 


dt 


— ln(l - zf) , 


t 


Sp{Q) = 0 , ^p(l) = ^ ^ , 

TT^ 

Sp{z) = -Sp{l - z) + — -ln( 2 :) ln(l - z) , 

D 

Sp(z) = -SpJ-)-fr-\\xf(-z). 
z b 1 

Useful identities for loops: 

arctan( 2 :) = — In , arctanh( 2 :) = - In ^ . 

2 ?- 1 %z 2 1 z 

Phase space element, dfp = \pf d\p\ d cos 9 df , cos 9 & [—,0e [0,27r[. 

= [ d'^p 6{p‘^ — TTi^) 9{E) ; E = s/fp -\-TTiJ . 

2E I 


Dirac algebra identities, for more see ref. [163], especially the appendix. 


2 \n1 JuJu 9fJ,U 


(A.23) 

(A.24) 

(A.25) 

(A.26) 

(A.27) 

(A.28) 

(A.29) 

(A.30) 

(A.31) 

(A.32) 

(A.33) 

(A.34) 

(A.35) 
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APPENDIX A. GENERALITIES 


A useful tool within this context is the TRACER routine [164] running under the symbolic alge¬ 
bra program mathematica. 

Chiral projectors L{R) = (1 =f 75 ) 72 : 

= 1 , 75 = 75 , (A.36) 

{L{R)f = L{R) , LR = RL = 0, {L{R))^ = L{R) . (A.37) 

Eurther we have 

72 = 1 , 7t = ^o. (A.38) 

Eermion fields : 

^ = (^*)^7o , = L{R)'il) , 

^ = L{^={L{R)ij)^-^, = i^\L{R))^^o = i^^loR{L)=^R{L) . (A.39) 




Appendix B 


The large energy expansion 


B.l Feynman Rules in the Large Energy Limit of QCD 


In the hadron limit of an infinitely heavy meson M with mass uim ^ oo and large energy for 
the final one Ep(^v) the QCD Lagrangian defined in (A.l) is just: 

= Cleet + C^(l/-E) , (B.l) 


where the effective Lagrangian is given by 

Cleet = 


- ■ n 

Qn — in--D q„ 


(B.2) 


where qn{x) = q{x) are the large components of the light quark spinor field. The 

Up = 2v — n_ is a light-like vector with = 2 and Eg ^ E is the energy of the light 

quark. Here v and n denote respectively the four-velocity = 1) of the heavy quark Q with 
momentum pq = mqv + k and the four light-like vector (n^ = 0) parallel to the the light quark 
Pq = En_ + k', where the small residual momentum k (and k' )of order Aqcd- 
It follows immediately, the Feynman rules of the LEFT formalism: 



n-k+i e 2 



where the light quark q is represented by a double line. 
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APPENDIX B. THE EARGE ENERGY EXPANSION 


B.2 The Factorizabe corrections 


We give here the 0{as )-factorizable corrections to the EEET form factors (defined in Eqs. (2.68)- 
(2.77)). We now present the result for the hard scattering correction to B ^ P form factors, as 
defined by A/_|_ o,r in Eqs. (2.68), (2.70). The renormalisation convention* implies A/+ = 0 by 
definition. Thus, they are then given by [63] 


A/+ = 

0 , 

(B.3) 

A/o = 

mB - 2Ep 

2Ep 

(B.4) 

A/p = 

mp + mp 

2Ep 

(B.5) 


with the quantity 


AFp = 


SF/s/p 1 

(/+ )+{u )p, 


Ncms 

is defined in terms of moment of the leading twist distribution amplitude (as usual u 


{u ^)p = du 


<P{u) 


u 


Eor the B meson, the quantity {l+)+ reads 


(B. 6 ) 

l-u) 

(B.7) 



(B.8) 


Note that A/o vanishes at s = 0 (Ep = mp/^) as required on general grounds. Eollowing the 
renormalisation convention, the hard correction to the other B ^ V form factors, defined in 
Eqs. (2.71)-(2.77), reads [63] 


AAi 

AA 2 


ATi 


AT 2 

AT 3 

We introduce the quantities 


AE = 0, 

my mp 


rriBirnB - 2 i?y) 


Ev tub - my 
rriB 


AEl 


AFii 


AEy 

- AF_l 

2 


AF, 


mp 


4F 


AF, + 


my / mp 


V 


Ey V 2F 


Vy 


AFi 


(B.9) 


(B.IO) 


*In ref. [63], it was convenient to define the factorisation scheme (or renormalisation conventions for the “soft 
form factors”) by imposing the condition that /+ = V = ^0 = exactly to all 

orders in perturbation theory. 
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where the two terms, namely {u ^)|| and {u ^)_l are given by 



/ 


du (l)\\{u)/u, 


and 



du 0 ±(m)/m. 


(B.13) 


(B.14) 


B.3 The Functions F- 

L 


In this appendix we list the functions and (defined in Eqs. (3.64) and (3.65)), 

representing the power correction to the matrix elements of the operators Og, Oi and O 2 respec¬ 
tively. The corrected one loop matrix elements of Og read [22]: 


and 



32 , u 8 s , . 8 . 4 11 - 16s+ 852 

-In-In s- ZTT -;-rr- 

9 mb 9 1-s 9 9 (1 - s)^ 

+ g _ ~^3 [(9s — 5s^ + 2s^) Bq(s) — (4 + 2s) Co(5)] ; (B.15) 



16 1 8 5-2s 

-Ins -|— -TTT 

9 1 - s 9 (1 - s)2 




where s = s/m^, Bo(S) = Bo(g^, ml) is given in (3.56), and the integral 


^(5) 



1 

a; (1 — s) + 1 


In 


X 


1 — x(l — x)s 


(B.17) 


can be expressed in terms of dilogarithms. The corresponding power correction to the matrix 
elements of the operators Oi and O 2 are given respectively (using = ln(/r/mb), Lg = ln(s) 
and rhc = me/mb ) as [22]: 


= 


1424 16 . 64 \ ^ 16 , , / 16 

+ 7777: + — Lc ] L^Lg + I 


729 243 


^ O ~ -4 

- — - Tfl 

2835 315 ^ 


32 


27 


LnS + 


f 16 




76545 


32 

-m, 

8505 


1215 135 

L~s^ 


■ m 


-2 


LnS 


^ 2 , .(9) 

243 F h , 


(B.18) 


and 


= 


256 32 . 128 ^ 

^ ~ ^ ~ “9~ 


32 

^+81 


32 ^ 64 , ~ 

■405 + 45’"- 


+ 


16 


945 105 


■ m 


-4 


s + 


32 


+ 


64 


25515 2835 


■ m 


-6 


L„ 
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APPENDIX B. THE EARGE ENERGY EXPANSION 


—v+/f 

81 ^ 


(B.19) 





(B.20) 


The analytic results for f['^\ f 2 ^\ and (expanded up to and {rhl)^) are rather 
lengthy. The formulas become relatively short, however, if we give the charm quark mass depen¬ 
dence in numerical form (for the characteristic values of ?fic=0.27, 0.29 and 0.31). 

We write the functions as 


Li (a = 1, 2; 6 = 7, 9; * = 0,..., 3; j = 0,1). (B.21) 

The numerical values for the quantities ka\i,j) are given in Tables (B.l) and (B.2). 




fhc = 0.27 


The = 0.29 


TTic = 0.31 


1 1 

1 

-12.327 

+ 

1 

0.13512i 

1 

-11.973 

+ 

1 

0.16371 i 

1 

-11.65 

+ 

1 

0.18223 i 


-0.080505 

— 

0.06718H 

-0.081271 

— 

0.059691 i 

-0.080959 

— 

0.051864i 

kf\l,0) 

-33.015 

— 

0.42492 i 

-28.432 

— 

0.25044 i 

-24.709 

— 

0.13474 i 


-0.041008 

+ 

0.0078685 i 

-0.040243 

+ 

0.016442 i 

-0.036585 

+ 

0.024753 i 

k?{2,0) 

-76.2 

— 

1.5067i 

-57.114 

— 

0.86486 i 

-43.588 

— 

0.4738 i 

kf\2,l) 

-0.042685 

+ 

0.015754 i 

-0.035191 

+ 

0.027909 i 

-0.021692 

+ 

0.036925 i 

M'^(3,0) 

-197.81 

— 

4.6389 i 

-128.8 

— 

2.5243 i 

- 86.22 

— 

1.3542 i 

M'^(3,1) 

-0.039021 

+ 

0.039384 i 

-0.017587 

+ 

0.050639 i 

0.013282 

+ 

0.052023 i 


Table B.l: Coefficients in the decomposition of for three values of fhc [22]. 


TTlc = 0.27 TUc = 0.29 TUc = 0.31 


4 ®'( 0 , 0 ) 

7.9938 

— 

0.81071 i 

6.6338 

— 

0.98225? 

5.4082 

— 

1.0934i 

4°'(o,i) 

0.48303 

+ 

0.40309 i 

0.48763 

+ 

0.35815? 

0.48576 

+ 

0.31119? 


5.1651 

+ 

2.5495 i 

3.3585 

+ 

1.5026? 

1.9061 

+ 

0.80843? 

4”(i.i) 

0.24605 

— 

0.047211 i 

0.24146 

— 

0.098649? 

0.21951 

— 

0.14852 i 

4 ®'( 2 , 0 ) 

-0.45653 

+ 

9.0402 i 

-1.1906 

+ 

5.1892 i 

-1.8286 

+ 

2.8428? 

4”(2,1) 

0.25611 

— 

0.094525? 

0.21115 

— 

0.16745? 

0.13015 

— 

0.22155? 

4 ®' (3.0) 

-25.981 

+ 

27.833? 

-17.12 

+ 

15.146? 

-12.113 

+ 

8.1251? 

fcf(3,l) 

0.23413 

— 

0.2363? 

0.10552 

— 

0.30383? 

-0.079692 

— 

0.31214i 

k?\0,0) 

4.3477 

+ 

0.56054 i 

4.0915 

+ 

0.44999? 

3.8367 

+ 

0.3531? 

kP{0,l) 


0 



0 



0 


kPifO) 

1.5694 

+ 

0.9005? 

1.4361 

+ 

0.73732 i 

1.3098 

+ 

0.60185? 

kPifl) 

0.0010623 

— 

0.12324i 

-0.016454 

— 

0.11806? 

-0.031936 

— 

0.10981? 

kP{2,0) 

-0.14311 

+ 

1.2188? 

0.011133 

+ 

1.05? 

0.13507 

+ 

0.89014 i 

kP{2,l) 

-0.12196 

— 

0.099636? 

-0.13718 

— 

0.068733? 

-0.14169 

— 

0.035553? 

kP{3,0) 

-2.5739 

+ 

0.5952H 

-1.6949 

+ 

0.76698? 

-1.0271 

+ 

0.77168? 

kP{3,l) 

-0.18904 

— 

0.0025554 i 

-0.17416 

+ 

0.049359? 

-0.13592 

+ 

0.093? 


Table B.2: Coefficients in the decomposition of f^'* and f 2 ^ for three values of fhc [22]. 





























Appendix C 


B K*i+i- in SUSY 


C.l The functions fi{x) 


In this appendix we give the various one-loop integral functions, namely fi{x) {i = 1, • • •, 8), 
which appear within the MSSM contributions to the magnetic dipole moment and to the semilep- 
tonic coefficients. They are given respectively by [135]: 


flix) 

—7 + bx + Sx^ 2x — 3a;^ , 

= 6(1 

(C.l) 

f2{x) 

3a; — 5a;^ 2x — 3a;^ 

2(l-a;)2 + (l-a;)3 

(C.2) 

fsix) 

2 + 5a; — a;^ x 

Q{l-xf (l-a;)4 

(C.3) 

U{x) 

1 + a; X 

2(l-a;)2 + (1-a;)3 

(C.4) 

f5{x) 

X X , 

= 1 

1 — a; (1 — xy 

(C.5) 

Mx) 

38x — 79x^ + 47a;^ 4a; — 6a;^ + 3a;"‘, 

6(1 

(C.6) 

frix) 

52 — 101a; + 43a;^ 6 — 9a; + 2a;^ 

6(1-:.)=> + (1-x)- 

(C.7) 

fs{x) 

2 — 7a; + lla;^ 6a;^ , 

^3 N4log^- 

(1 — a;)'* (1 — xy 

(C.8) 
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C.2 The Auxiliary Functions Ci{m\^ m^), di{m\^ ml) 


The various auxiliary functions Ci{m\, ml, ml) and di{m\, ml, ml, m\) are listed below [135]: 


/ 2 2 2 \ 
Co(mi,m2,m3) 

C2{ml,ml,ml) 


2 i 

mf log ^ 




L (ml — ml){m{ — ml) 


2 1 "^1 
mf log — 


+ (mi <->• m 2 ) + (mi m 3 ) 


, (C.9) 




(mf — ml){m\ — ml) 


+ (mi <->• m2) + (mi m3) 


(C.IO) 


do{ml,ml,ml,ml) = - 


2 i "^1 

mi log — 

__ 

L(mi — m 2 )(mi — m|)(mi — m|) 


+ (mi ^ m 2 ) 


+ (mi <->• m3) + (mi <->• mi) 


(C.ll) 


^ / 2 2 2 2 \ 
d2[m^,m2,m^,m^) = — 


4 1 

m\ log ^ 


4 L(mf — ml){m\ — ml){m\ — ml) 


+ (mi ^ m 2 ) 


+ (mi <->• m 3 ) + (mi mi) 


(C.12) 
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